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Abstract 

Both Grand Unified symmetries and discrete flavour symmetries are appealing 
ways to describe apparent structures in the gauge and flavour sectors of the Standard 
Model. Both symmetries put constraints on the high energy behaviour of the theory. 
This can give rise to unexpected interplay when building models that possess both 
symmetries. We investigate on the possibility to combine a Pati-Salam model with 
the discrete flavour symmetry S4 that gives rise to quark-lepton complementarity. 
Under appropriate assumptions at the GUT scale, the model reproduces fermion 
masses and mixings both in the quark and in the lepton sectors. We show that in 
particular the Higgs sector and the running Yukawa couplings are strongly affected 
by the combined constraints of the Grand Unified and family symmetries. This in 
turn reduces the phenomenologically viable parameter space, with high energy mass 
scales confined to a small region and some parameters in the neutrino sector slightly 
unnatural. In the allowed regions, we can reproduce the quark masses and the CKM 
matrix. In the lepton sector, we reproduce the charged lepton masses, including 
bottom-tau unification and the Georgi-Jarlskog relation as well as the two known 
angles of the PMNS matrix. The neutrino mass spectrum can present a normal or an 
inverse hierarchy, and only allowing the neutrino parameters to spread into a range 
of values between A -2 and A 2 , with A ~ 0.2. Finally, our model suggests that the 
reactor mixing angle is close to its current experimental bound. 
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1 Introduction 



Gauge coupling unification suggests that the Standard Model (SM) gauge group is gener- 
ated when a unified larger gauge group is broken at a very high energy scale compared to 
the electroweak (EW) one. In this picture SM fermions are accommodated in representa- 
tions of the unified gauge group G and an appropriate scalar Higgs sector is introduced both 
to trigger the spontaneous breaking of G down to SU(3)c X U(l) em and to reproduce the 
fermion mass matrices. Typically Grand Unified Theories (GUTs) are quite constrained 
and non-trivial relations are obtained among the SM fermion mass matrices. On the other 
hand it is a common belief that fermion mass and mixing matrix patterns can be explained 
by appealing to a flavour symmetry Gf and in the last years particular attention has been 
devoted to the study of discrete flavour symmetries, thanks to their simplicity in recovering 
realistic lepton mixing patterns. 
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Table 1: Neutrino oscillation parameters from two independent global fits 

The present neutrino oscillation data [l}j3] are summarised in table [TJ where we display 
the results of two independent global fits. The pattern of the mixings is characterised by 
two large angles and a small one: 8 l 23 is compatible with a maximal value, but the accuracy 
admits relatively large deviations; 9[ 2 is large, but about 5a far from the maximal value; 9[ 3 
has only an upper bound. According to the type of the experiments which measured them, 
the mixing angle 8 l 23 is called atmospheric, 9\ 2 solar and 8[ 3 reactor. We underline that there 
is a tension among the two global fits presented in table [T]on the central value of the reactor 
angle: in jlj we can read a suggestion for a positive value of sin 2 9[ 3 ~ 0.016 ± 0.010 [1.6<r], 
while in pj a best fit value consistent with zero within less than la is found. Therefore 
we need a direct measurement by the future experiments like DOUBLE CHOOZ |4|, Daya 
Bay D, MINOS §, RENO |7|,T2K § and NOvA @. 

The closeness of the leptonic atmospheric angle 9 23 to the maximal value gives relevant 



indications on the flavour symmetry: it is well known 10[ 11 that a maximal 8 l 23 is not 



achievable with an exact realistic symmetry. This forces to study models based on the 
breaking of the flavour symmetry and a promising choice is the kind of realizations based 
on non-Abelian discrete groups which reproduce the lepton mixing matrix of the so-called 
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tribimaximal (TB) pattern [12) (sin 2 0f 2 B = 1/3, sin 2 6^ = 1/2 and sin6>f 3 B = 0) at lead- 
ing order. This mixing scheme represents a very good approximation of the experimental 
measurements: the TB values for the atmospheric and the reactor angles are inside the 
la error level, while that one for the solar angle is very close to the upper la value. The 
corrections from the symmetry breaking provide perturbations to the angles and in partic- 
ular a deviation from zero for the reactor angle, in agreement with the recent indication 
of a positive value for 9[ 3 [l]. Despite this success of the TB mixing scheme, there are 
some complications, in particular when combining the flavour symmetry with GUTs. In 
the following, we will comment on this and motivate the viability of a competing scheme, 



the bimaximal (BM) mixing pattern 13 



A lot of effort has been put in reproducing the TB pattern by the use of non-Abelian 
discrete symmetries: the best known groups implemented in the construction of flavour 



models are A4 14,15 , S4 [16] , T' 17 and A(27) [18] . A common feature among many of 
these realizations is to get a spontaneous breaking scheme responsible for the TB mixing 
by the use of a convenient assignment of the quantum numbers to the SM particles and the 
introduction of a suitable set of scalar fields, the "flavons" , which, getting non-zero vacuum 
expectation values (VEVs), are responsible for the symmetry breaking of Gf. A central 
aspect of the model building is the symmetry breaking chain: the flavour group is broken 
down to two distinct subgroups, which correspond to the low-energy flavour symmetries of 
the charged leptons and of neutrinos. 

Further studies have been presented in which the SM is extended to a GUT scenario: 
in our opinion, these analyses emphasize the difficulty in the construction of a flavour 



GUT model 19 which naturally lead to realistic fermion phenomenology and to a fair 



explanation of the gauge symmetry breaking chain to get SU(3)c x U(l) em . Concerning 
the case of the flavour GUT models 20 that predicts the TB mixing in the lepton sector, 
the combinations of the constraints arising by the flavour symmetry and by the GUT 
group lead to wrong predictions for the fermion masses and mixings. The problem is 
usually avoided by recurring to non-minimal Higgs or flavon field content and by assuming 
peculiar symmetry breaking patterns of the GUT gauge symmetry and ascribing quite 
often at type-II See-Saw as the origin of the neutrino masses. Moreover these patterns are 
often not supported by the study of the scalar Higgs potential, leaving open the question 
if such peculiar patterns may be actually realized or not. 

Moreover, a common feature of flavour models which deal with the TB scheme is a 
value for the reactor angle very close to zero, in the absence of specific dynamical tricks 
(see (21 for a model in which such a trick is implemented). However, if the next future 



neutrino-appearance experiments will find a value for 8[ 3 close to its present upper bound, 
about the sine of the Cabibbo angle A, the TB mixing should be considered an accidental 
symmetry. In this case a new leading principle would be necessary, for which quark-lepton 
complementarity [22], 9[ 2 + A ~ 7r/4, would be a good candidate. This has the advantages 
that it can be naturally related to GUTs and does not necessarily predict a small reactor 
mixing angle. 



This idea was developed in a model 23 based on the S4 discrete group. In this con- 
struction the PMNS matrix coincides with the bimaximal mixing 13 (sin 2 #^ A/ = 1/2, 



sin 2 02 3 M = 1/2 and sin9f 3 M = 0) in first approximation; since the BM value of the solar 
angle exceeds the 3a error, large corrections are needed to make the model agree with 
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the data; the perturbations are naturally constrained to get the "weak" complementarity 
relation, 6[ 2 + 0(X) ~ 7r/4, and sin#^ 3 ~ A in most of the parameter space. The model 
only deals with the lepton sector and an extension in the quarks sector is lacking. In this 



paper we aim at revisiting the model in 23 in order to include a realistic description of 
quarks. 

Our starting point lies in the complementarity relations: 

6[ 2 + 6\ 2 ~ tt/4 + 0(X 2 ) , 4 + 9 q 23 ~ -vr/4 + 0(X 2 ) . (1) 

For the third mixing angles we know that in the quark sector it is very small, 8f 3 = 0(A 3 ), 
while in the lepton sector, as already mentioned, it has only an upper bound, sin6^ 3 < A. 
We will see that in the model described below, we predict it to be 0(X). 

In non-GUT contexts, no compelling model leading to the exact complementarity has 
been produced so far and indeed in |23| a weaker version of eq. ([I]) has been used in which 
the quark contributions are substituted by similar terms originating from the charged 
lepton sector. On the other hand, exact complementarity is possible in cases where the 
flavour symmetry group is combined with a GUT group, as in the Pati-SalamQ context 24 , 
where the charged lepton and the down-quark mass matrices are similar, 

M e ~ M d . (2) 

Other popular GUTs are SU(5) and 50(10), which however are less appealing when 
trying to recover the QL complementarity relations. Indeed, in the minimal SU(5) [25] one 
has M e ~ Mj and as result, a correction of order A to the solar angle does not correspond 
to the Cabibbo angle of the CKM matrix. 

On the other hand a reason to prefer Pati-Salam over SO (10) is related to the type-I 
and type-II See-Saw mechanisms. In these two GUT contexts, we expect both left-handed 
(LH) and right-handed (RH) neutrino Majorana mass matrices to be present. As a result, 
the effective LH neutrino mass matrix will get the contributions through the type-I as well 
as the type-II See-Saw mechanisms. In general, and this happens also in our proposal, 
this interplay introduces two mass scales and a highly non-trivial flavour structure for the 
effective neutrino mass matrix, which difficultly leads to a realistic description of the PMNS 
matrix. For this reason, a hierarchy between the two contributions is usually assumed. As 
we will see we can reproduce quark-lepton complementarity in our model, if the type-II 
See-Saw is dominant. 

This possibility has already been investigated in several flavour GUT models, for exam- 



ple in 26 in the context of the 50(10) GUT. However, in 27 it has been agued that the 
type-II dominance in the context of minimal 50(10) models is highly disfavoured. Even 
so, when restricting to particular supersymmetric parameter space, the type-II dominance 
could be possible both in the minimal and non-minimal 50(10) approaches, [28). In the 
Pati-Salam context, there is much more freedom and the eventual dominance of one of the 
two contributions could be realized. In this paper we study the gauge Higgs potential and 
we verify that a type-II dominance can be justified, even if it puts strong constraints on 
the model building realization. 

Adopting a Pati-Salam context, we loose the exact gauge coupling unification, but at the same time 
we avoid several complications which are instead present in other GUT scenarios, as described in the text. 
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After a detailed analysis of the fermion phenomenology, we move to the study on the 
Higgs gauge sector, which is responsible for the gauge symmetry breaking steps to finally 
get SU(3)c x U(l) em . Thus we consider the renormalization group equations (RGEs) both 
for the gauge couplings and for the fermion masses and mixings from the cutoff of the 
theory down to the low-energy scale. Here we anticipate that the RGEs analysis is crucial: 
the model turns out to be viable only in a small region of the parameter space. On one 
side the gauge coupling RGEs analysis is deeply modified by the non-minimal Higgs field 
content required by the presence of the flavour symmetry Gf and puts strong constraints 
on the scales of the model. On the other side the GUT flavour mass matrix structures 
interfere with each other, through the Yukawa RGEs, with non-negligible consequences for 
the neutrino phenomenology. 

The paper is organized as follows. In section [2] we describe the flavour structure of 
the fermion mass matrices in order to recover realistic lepton and quark mixing matrices 
and fermion mass hierarchies. In section [3] we enter in the details of the model building 
construction, specifying the transformations of all the fields of the model under the gauge 
and flavour groups and discussing the flavon vacuum misalignment necessary to the flavour 
symmetry breaking chain. Afterwards we deal in section [4] with the study of the gauge 
Higgs potential and we analyze the constraints coming from the flavour symmetry on the 
scalar field content. In section [5] we perform the analysis of the renormalization group 
running of fermion masses and mixings and of the gauge coupling constants from the GUT 
scale down to the low-energy scale. In section [6] we report the phenomenological analysis 
after the running evolution. Finally in section [7] we conclude. In the appendices we report 
technical details. 



2 Outline of the model 

In this section we present a general discussion on the flavour structure of fermion masses 
and mixings resulting form the analysis of the complementarity relations. Eq. ([I]) suggests 
that the angles in the CKM and PMNS matrices may have a common origin. Looking at 
their definitions 

V = VtV d , U = U\U U , (3) 

where V u , Vd, U e , U u diagonalize M u Ml, \ /,/ \ V d . M e M\ and m v respectively, we see that 
this common origin can be motivated in a GUT context where some relations among the 
mass matrices are present. In PS models, we have the following expression which links 
down-quarks and charged leptons 

U e ~V d ; (4) 
we will see in a while how this enters in the model construction. 



In 23 the two large lepton angles in the PMNS matrix arise only from the neutrino 
sector: in the first order approximation, both the solar and the atmospheric angles are 
maximal and after that some corrections from the charged lepton sectors lower the value 
of the solar angle in order to accommodate the data. In our realization one maximal angle 
originates from the charged leptons and the other from the neutrinos: afterwards some 
corrections are necessary to make the solar angle agree with the measurements and we 
require that these corrections come from the charged leptons. This is an other possible 
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choice with respect to 23 and in the next section we will show that it is easily realized in 
our context. 

In the approximation of small A, the lepton mixing matrix can schematically be written 



as 



U = # 23 (-|) #i 3 (A)#i 2 Q - A) = (# 23 (-|) #i 3 (A)# 12 (-A)) # 12 Q , (5) 



u 



where Rij(0) stands for a rotation in the (ij) plane of an angle 6. ^ As a result, after a 
suitable commutation of matrices, U e can naively be written as 

U e = #23 (|) #13 (A) #12 (A) . (6) 

The CKM matrix is given in first approximation as 

V = R 12 (X). (7) 

Because of eq. Q, Vd has the same structure as U e in eq. ^ and therefore we can obtain 
the CKM matrix as 

V = (#i2(-aA) J R 13 (-A) J R 23 (-^)) J R 2 3(^) J Ri3(A) J Ri2(A)) . (8) 



^ .. ' ^~ 



We see that the angles of the (23) and (13) rotations in V£ should be the opposite of 
those in Vd, while the angles in the (12) sector should be different. We have schematically 
indicated this via the a coefficient. Analogous to eq. ([6]), we write V u as 

V u = #23 Q) #i 3 (A)#i 2 («A) . (9) 

Moving to the explicit form of the mass matrices, the generic effective Majorana neu- 
trino mass matrix m u which is diagonalized by U u as in eq. (§, U„ = #12 (f), through 
m diag _ jjj m^JJu, is given by 



m v ~ lb a . (10) 




In the charged lepton sector, we are looking for the most general form of the mass 
matrix, whose square M e M\ is diagonalized by the action of U e of eq. (|6]), 

{M e Mt) dia9 = UlM e M\U e . (11) 



Inverting eq. (11 ) we find in the limit m e — > 



2 

mi 





A 






j 


1 


i 




1 





M e Ml~-±\ A 1 1 + — ^ A 1 -1+..., (12) 




2 A coefficient of order one is in general present to multiply each angle, but we do not show them here 
in order to simply the discussion. We will consider the precise expressions in the following sections. 
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where the dots stand for suppressed contributions. This can be obtained if M e is given by 




M e p 1 + 1 +.... (13) 



It is interesting to note that, moving to the basis of diagonal charged leptons and consid- 
ering only the leading order terms, the neutrino mass matrix results to be of the BM type, 
already presented in the introduction. 

As already stated, the PS gauge structure implies the relation M e ~ Md and therefore 
the down-quark matrix has a similar structure as in eq. (13). Looking at eq. ^ we find 



that also M u should have a similar structure, but here there is a slight difference due to 
the a coefficient. 

In order to reproduce the result in eq. ([8]), i.e. exact cancellations in the (23) and (13) 
sectors but not in the (12) sector, we need to study the origin of the different rotations in 
Q. It is easy to see that in the diagonalization of the mass matrices in eq. (12), the (12) 



rotation originates from the second families, while the (13) rotation comes from the third 



families. We conclude that the up-quark mass matrix has the same form as (13), in which 
the third column is proportional to that of the down-quarks, while the second column is 
not. 

As already stated in the introduction, we expect contributions on the effective neutrino 
mass matrix coming from both the type-I and the type-II See-Saw mechanisms. For the 



moment we just assume that the type-II is the only responsible of eq. (10) and we verify 
in section [4] that it is indeed dominating with respect to the type-I terms. 

In the next section we enter in the details of the model building, explaining the origin 
of the mass matrices displayed above. 



3 The flavour model building 

In this section we define our framework. The model is based on the PS gauge symmetry, 
SU(4)o x SU{2) L x SU(2)r, present at high energy where a supersymmetric context is 
assumed. The complete flavour group Gf is the same as in 23 given by the product of 
the following different terms: 

G f = S 4 xZ 4 x U(1) FN x U(1) R . (14) 

The group 5*4 is the permutation group of four distinct objects, isomorphic to the 
group O which is the symmetry group of a regular octahedron. It has 24 distinct elements 
filled in five conjugate classes and therefore it has five irreducible representations, two 
one-dimensional denoted as li and I2, one two-dimensional labelled as 2 and two three- 
dimensional written as 3i and 32- Here we recall the multiplication rules, while the Clebsch- 
Gordan coefficients, the explicit structures of the generators, the list of the elements and 
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further details on the group theory of S4, are reported in appendix [A| 



11 x R — R x li — R where R stands for any representation 

1 2 x 1 2 = li 
1 2 x 2 = 2 
1 2 x 3i = 3 2 
1 2 x 3 2 = 3i 



2x2 = l! + l 2 H 
2 x 3i = 3i + 3 2 
2 x 3 2 = 3i + 3 2 



(15) 



3 1 x 3 1 = 3 2 x 3 2 = l 1 + 2 + 3 1 + 3 2 
3i x 3 2 = 12 + 2 + 3! + 3 2 . 

The spontaneous symmetry breaking of S4 is responsible for the flavour structure of 



the matrices in eqs. (10, 13): S4 is broken down to two distinct subgroups and it is the 



presence of this mismatch at the LO which allows to construct eqs. (10, 13). More in 
detail, the different subgroups to which S4 is broken down are the subgroups preserved by 
the VEVs of the flavons. In order to determine these structures, it is necessary to identify 
the elements of the group which leave the VEVs of the flavons invariant under their action. 
Doing so, we find that 5*4 is broken down to Z 2 x Z 2 in the neutrino sector, originated by 
the elements ST 2 S and T 3 ST of the classes C 2 and C3, respectively. In the charged fermion 
sector, we anticipate that the superpotential is invariant under an accidental Z 2 symmetry 
and therefore the identification of the residual subgroup must account of this additional 
term. We find that S4 is broken down to a Z 2 x Z 2 group, distinct to that one in the 
neutrino sector, generated by the two elements — T 2 and —TST 2 S of the classes C 2 and 
C 3 , respectively. One can argue that the group Z 2 x Z 2 implies a degeneracy among two of 
the charged fermion families, but this does not happen, because the two Z 2 relate different 
families: the accidental Z 2 exchanges F 2 and while the Z 2 subgroup of S4 relates F± 
and F 2 . 

The other terms in Gf carry out other roles: the Abelian Z 4 symmetry is required 
to avoid dangerous terms in the superpotential, it helps to keep the different sectors of 
the model separated, quarks from leptons and neutrinos from charged leptons, and it is 
also useful to guarantee the flavon vacuum alignment; the continuous Frogatt-Nielsen (FN) 
Abelian symmetry [29], U(1)fn, is introduced to justify charged fermion mass hierarchies; 
the continuous i?-symmetry U(1)r, that contains the usual .R-parity as a subgroup, is a 
common feature of supersymmetric formulations and simplifies the constructions of the 
scalar potential. It is worth stressing that the supersymmetric context is of great utility in 
the discussion of the scalar potential and helps in the gauge coupling running, but a similar 
non-supersymmetric model can be constructed as well. In particular in the present paper 
we only deal with SM particle and therefore we use the same symbols for a supermultiplet 
and its even .R-parity components. 
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3.1 The matter, Higgs and flavon content of the model 

In the PS context, the five matter multiplets of each family of the SM plus a RH neutrino 
and their superpartners are unified in only two supermultiplets: a LH and a RH ones as 
follows, 

SU(3) C x SU(2) L x U(1) Y -> PS 

(3, 2, 1/6) Q + (1, 2, -1) L -> (4,2,1) (16) 

(3, 1, -2/3)„c + (3, 1, l/3) d c + (1, 1, l) e c + (1, 1, 0) v c (4, 1, 2) 

The three copies of the LH supermultiplet are combined in the three-dimensional repre- 
sentation 3i of S4, while the three families of the RH supermultiplet are in 1 2 , I2 and 
li respectively. The fact that we can put different representations within one family in 
different representations of the family symmetry group is essential here. Note that this 
would not be possible in (minimal) SO (10), where all Standard Model particles are in one 
sixteen dimensional representation. The first two families are also charged under U(1)fn 
by two units. This suppresses their masses with respect to the third family ones. Further 
suppression of the first family with respect to the second is due to their different Z4 charges. 
All the properties of the matter fields are summarized in table [2] 



Matter 


F L 








PS 


(4,2,1) 


(4,1,2) 


(4,1,2) 


(4,1,2) 


S, 


3i 


1 2 


1 2 


li 


U(1) FN 





2 


2 





Z4 


1 


1 


i 


—i 



Table 2: Transformation properties of the matter fields. Notice that the PS assignments should be read in 
agreement with SU(4) C x SU(2) L x SU{2) R . 

Our model contains five flavon fields: two S4 triplets (<p and ip') that, because of their 
Z4 charge, deal at LO only with the Dirac Yukawa couplings of quarks and leptons, and two 
fields, one singlet (a) and one triplet (x), that, by Z4, deal at LO only with the Majorana 
masses of neutrinos. The fifth flavon is the Froggatt-Nielsen messenger, which we indicate 
with 9. Their properties are shown in table [3j Under the continuous .R-symmetry, the 
matter superfields transform as U(1)r = 1, while all the flavons are neutral. 

Fermion masses and mixings arise from the spontaneous breaking of the flavour symme- 
try by means of the flavons which develop VEVs according to the following configuration: 
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Flavons 


e 


<p 


<p' 


X 


<7 




ii 


3i 


3 2 


3i 


ll 


U(1) FN 


-i 














z 4 


1 


% 


% 


1 


1 



Table 3: The flavon field content and their transformation properties under the flavour symmetries. All 
fiavon fields are singlet of the gauge group. 



at LO we have 

(<p) = I i K . (^') = I 1 I v > ( 17 ) 





(X> = k, (a) 



(0) = v B . (19) 
In this section we simply assume this VEV alignment and we will prove it to be a natural 



solution of the minimization of the scalar potential in section |3.4| Furthermore we assume 
that the FN messenger and the other flavons have VEVs of the same order of magnitude: it 
results partly from the minimization procedure and partly from the constraints coming form 
the comparison with the measured mass hierarchies, as it will be clearer in the following. 
Without loss of generality it is useful to keep the notation compact and write 

VEV 

« A , (20) 

where VEV refers to the vacuum expectation value of any flavon of the model and A ~ 0.2 
is close to the Cabibbo angle. A is the flavon cut-off that we assume to be the largest scale 
in the model: it corresponds to the scale of the flavon dynamics. 

The insertion of flavons in the mass terms, leads to non-renormalizable operators, as 
we will see in the following section. This means that we can write the superpotential as an 
expansion in powers of flavon/ A and we can stop the expansion after the first orders. 

The Higgs fields of our model relevant to build the fermion mass matrices transform 
under the gauge group and under the Z 4 factor of the flavour symmetry: in table [4] we 
summarize their transformation rules. The first three fields, 0, and p, deal at LO only 
with the Dirac Yukawas. Due to the Z4 charges, <fi and 0' are responsible of the third family 
and the charm quark masses, while p is responsible for the strange and p masses. The field 
p ~ (15, 2, 2) being in the adjoint of SU (4)c may develop VEV along the SU (4)c direction 
diag(— 3, 1, 1, 1). This implies that the leptons which get mass via this field are a factor 
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3 heavier than the corresponding quarks and therefore this field is very useful to describe 
the second family, at least in the down sector, reproducing the well known Georgi-Jarlskog 
relation [30], m 3m s , at the high energy scale. As we will see in the next sections, 
in order to recover the up-quark mass hierarchies the p projection along the light doublet 
Higgses, Dp and v^, has to point only in the down direction: the requirement = can 
be realized only if the Higgs field content contains two identical copy of the Higgs field 
(1, 2, 2) and this justify the presence of and 0'. 

Finally, as we will see in detail in the following sections, the field Ar is necessary to 
conclude the PS symmetry breaking pattern and to recover the SM gauge group through 
its spontaneous symmetry breaking VEV. At the same time, when Ar develops a VEV, 
it gives a Majorana mass to the right-handed neutrinos thus contributing to the effective 
neutrino mass matrix through the usual type-I See-Saw mechanism. A second source for 
the neutrino mass matrix comes from Al, in terms of the type-II See-Saw mechanism. 



Higgses 


0,0' 


P 


A L 




PS 


(1,2,2) 


(15,2,2) 


(TO, 3,1) 


(10,1,3) 


z 4 


1 


-1 


1 


-1 



Table 4: The Higgs fields responsible of generating fermion mass matrices and their transformation under 
the gauge and the Abelian flavour symmetries. All Higgs fields are singlets under S4 x U(1)fn X U(1)r, 
while they can transform under the Z4 factor. 



In our scheme the neutrino mass matrix is dominated by type-II See-Saw. As already 
stated in the introduction, the flavour structure of the effective neutrino mass matrix in 
eq. (10) arises from an interplay between the two See-Saw sources. In the PS context 
vtld ~ M u and this suggests a hierarchical structure for the type-I contribution, which does 
not agree with the flavour structure in eq. (10). We will show that, given rrin ~ M u , 
the required flavour structure for the Majorana RH neutrino mass necessary to recover eq. 
( 10 ) is not allowed in our model. This suggests to find a construction in which the type-II 
contributions are dominating over the type-I and we will show that such a feature puts 
strong constraints on the model building. 



3.2 Fermion mass matrices at leading order 

In this section we study the matter superpotential, Wy, and the resulting mass matrices 
for all the fermions. The Yukawa superpotential can be written as a sum of two pieces, 

W Y = W Dir + W Maj (21) 

containing respectively terms with Dirac and Majorana structures. Furthermore, making 
a power expansion in terms of flavon/A, we distinguish between leading and subleading 
couplings. We refer as "leading order" operators to those terms which provide M e ^, u and 
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m v as in eqs. (10, 13) in the limit A — > 0. The subleading orders corresponds to operators 
of higher dimensions. In this way it is easier to identify each entry of the mass matrices 
with the corresponding operator in the superpotential and furthermore it underlines the 
relevance of the subleading contributions. 

3.2.1 Dirac mass terms 

We first study the Dirac matter superpotential at LO which reads 

W L D ° = ift^W + <P)<P+ 

1 4 1 

+ y 2 -F L F°e 2 p( p' + y mT ,F L F^{ct> + 0')xf } + (22) 

i=l 



1 1 

+ E vm^^^ + y^ F ^ 2 (d> + 0V 



i=i 



(1 21 

Here we use a compact notation to avoid the proliferation of coefficients: the term X\ ' 
represents a list of products defined as 

(23) 

X\ 2) = {<?, 

where each term represents all the different S4 contractions which can be constructed with 
those flavons; furthermore we indicate with y%((j) + 0') the combination <f> + nf V and 
similarly for y Zt{i) and y 5 . 

When the flavour symmetry is broken the model describes a non-minimal PS model in 
which the Yukawa couplings present a well defined structure. Then, when the PS gauge 
symmetry is broken to the SM gauge group 0, 0' and the colour singlet component of p mix 
in four light Higgses, two up-type and two down-type, h u ^ and h' u d . Thus at the EW scale 
0, 0' and p have non-vanishing projections to the light Higgs components that acquire a 
VEV. We will indicate these components as v^' d , v^f and v u p ' d . As already said, we need 
to impose that the p field has no projection along the up direction: v u p = 0. This can be 
realized because in terms of the light Higgs up-VEVs, i>" = (h u ) and v% = (h' u ), v u p is given 
by 

v u p = U 13 v^ + U 23 v^ (24) 



where the matrix U is introduced in the appendix [Bj The constraint = —U13/U23V™ can 
be imposed thanks to the freedom we have in the superpotential and in the soft potential. 



n Note that we could relax the condition v p — allowing a mild hierarchy between v p 



3 This requirement, imposed by hand, could be motivated by some symmetry argument, but to introduce 
a mechanism that could explain v p = 0, it would be necessary deeply modifying the structure of our model. 



In the present paper we just assume this fine-tuning and we refer to the Appendix B.2.1 for further details. 
We just anticipate that the fine-tuning we introduce in the model is similar to the fine-tuning which is 
universally accepted whenever the MSSM has to be recovered. 
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and Vp, for example of order A 2 , without affecting the final mass hierarchies, but in the 
following, for simplicity, we work under the assumption that = 0. 
The final Dirac fermion mass matrices we get are given by 



M e LO = -3 y 2 v d p \ 3 + y x , (25) 







M L d ° = [ y 2 ]v d p X 3 + [ y x J c[[\ . (20) 

. 2/i 

/ \ 

M u i0 = m£° = y 3 t^A 5 + I yi U«A , (27) 












( 







! 


2/3 





^A 5 + | 








;, 




-2/3 






Vo 





2/i / 



where we used the compact notation yivt to indicate y^ vtj d + yf > v^, d and absorbed 
all the non-relevant CG coefficients. Note that y% is the sum of the different yz,{i) and 
that, by construction, 7/1,2,3 can be considered complex coefficients with modulus of order 
1. Note also that the different numerical coefficients between charged leptons and down- 
quarks originate from to the presence of p instead of <fi (</>') in the superpotential. The 
operators which should give contributions to the first families (those proportional to y& 
and 7/5) are vanishing, thanks to the special flavon VEV alignment. As a final comment, 
we are neglecting at this level of approximation the contributions to ® coming from the 
operators proportional to y^: these terms, which preserve the anti-alignment of the second 
and third entries of the second columns, are A 2 suppressed with respect to the LO ones 
proportional to y 2 . 



In order to identify the mass matrices in eqs. (25)-(26)-(27) with those in eqs. (10 



13), we need to define the (complex) fermion masses as follows: 

m M = -3y 2 Vp\ 3 , m T = yiV^X , 

m s = y 2 v d X 3 , m b = y x v JA , ( 2 8) 

m c = y?,v^X b , m t = yiv^X . 



Now we can comment on the masses as well as on the mixing matrices which can be driven 
at this approximation level. Notice that the top-quark Yukawa does not come from a 
renormalizable coupling and it presents the same suppression as the other third family 
fermion masses: as a result the dominance of the top-quark mass can be justified by the 
hierarchy between and vi. The mass matrices in eqs. (25)-(26)-(27) are diagonalized 
by a maximal rotation in the sector (23), i.e. U e = Vd = V u = R 
mass hierarchies are given by 



23 1 



vr/4) 



while the fermion 

















m b 



0(X 2 



m t 



0{X 4 



(29) 
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Furthermore, at the cutoff, we recover some relations among the masses of different fermions: 
the b — t unification and the Georgi-Jarlskog relation 30 



\m b \ 



\m. 



3|m s 



(30) 



Finally we should comment of the relative value of the top and of the bottom masses: 



m t 




(1) u , ( 2 ) u 

y \ •»% + y \ v p 






y { i )v i + y { i )v i' 



(31) 



Notice that usally this ratio is proportional to tan/3, the ratio between the up- and down- 
VEVs of the light Higgses, but this is not the case: indeed the light Higgses are combina- 
tions of 0, cf)' and p and therefore we may define tan /3 as 



tan (3 = 



y/(^) 2 + (v;,y 
K) 2 + K) 2 + W 



m t 



(32) 



3.2.2 Majorana mass terms 

We now discuss the part of the superpotential which contains the Majorana couplings. At 
LO0 it is given by 

= h F L F L A L + £ F L F L A L X® + ]T ^i^A.xf + z x F<F*A R , (33) 

i=l i=\ 

where we used the compact notation 

X (3) = {x,^}, 

(34) 

This superpotential is responsible for giving the following Majorana LH and RH neutrino 
mass matrices: 

/ k k t X \ / \ 

M L = I h\ k \v L , M R = \ v R , (35) 

\ k + k 2 \ 2 J \ z x ) 

where ^0,1,2 and z\ are coefficients of order 1 and vl, vr are the VEVs of respectively. 
In particular fc ,i,2 are defined in terms of k , k\(i) and £>2,(i): 

k = k + fci,( 2 )A + A; 2 ,(i)A 2 + & 2 ,(3)A 2 , 

h = fei,(i) + £2,(2) A , (36) 
fc 2 = -k 2 ,(i) ■ 

4 Regarding the terms which contribute to Mr, we consider at the LO only the first non vanishing term 
in the superpotential. 
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While Ml corresponds to the type-II contribution to the effective neutrino mass matrix, 
Mr provides a type-I term. Even at this approximation level, we note the tension between 
the two See-Saw contributions: m£ type " II ' ) = Ml presents a maximal rotation in the (12) 
sector, while it is easy to verify that ml typc_I ^ = mnM^m^ shows a democratic structure 
on the (23) sector which corresponds to a maximal rotation in this sector. To recover 



the mass matrix in eq. (10) it is necessary that the type-II contribution dominates over 
the type-I terms: in this case it is sufficient to identify a with ko, b with ki\ and c with 
k + k 2 \ 2 . 



The neutrino masses can be written as 



1 |2 

|mi, 2 | = 


(N 2 


T2\k \ 


\k\\ cos( 


9k 


|"7, 3 | = 


(m 2 


+ 2\k \ 


\k 2 \ cos( 


0k o 



OkJX + lhfX^vl, (37) 
e k2 )\ 2 )v 2 L , (38) 



with 9 ki the argument of the complex number fcj. We need cos(#fc — 9^) > in order 
to have \mi\ smaller than \ui2\- We see that in most of parameter space the spectrum is 
quasi degenerate, as the term with |A; | 2 that appears in all three masses dominates over 
the other terms that are A or A 2 suppressed. 

We also see that in most of the parameter space \m^\ is the central eigenvalue: indeed 
\ m \\ {\ m 2\) is shifted down (up) from the central value |fco| 2 f| by a term proportional to A, 
while |m 3 | stays closer to this central value as it is only shifted by a term proportional to 
A 2 . Having 1 777,3 1 as the central eigenvalue is obviously in contradiction with experimental 
data. We conclude that our model is only viable when the term 2\ko\ \k\\ cos (6^ — 9k ± )X is 
suppressed. This is clearly possible if either |fc | or \k\\ or cos(^ — 9 kl ) (or a combination of 
them) is small. In particular, the latter condition means that ko and k\ are almost perpen- 
dicular in the complex plane. We now investigate on these different scenarios, calculating 
the solar and atmospheric mass squared differences: taking as definition of Am 2 im the 
mass squared difference between the heaviest and the lightest neutrinos, we have different 
results for normal (NO) and inverse (10) mass ordering as given by 

Am 2 oi = m 2 , — m 2 



4|OA;i|Acos(0 feo - 9 kl )v 2 



k - Vkx)V L > 

2 



Am ltm 



Am atmNO = ( m 3 ~ m l) 
Am ltrnIO = ( m 2 - m l) 

2|fc o |Ncos(0 fco -9 kl )Xv 2 L T(\k 1 \ 2 - 2|/cop 2 |cos(0 Ao -9 k2 ))\ 2 v 2 L , 



(39) 



Am L/ 2 

On the right-hand side of the last equation, the first term is suppressed by A, while the 
second one by A 2 , and therefore we could conclude that it is the dominant contribution: it 
is however exactly the term which must be suppressed in order to avoid |rrz 3 1 as the central 
eigenvalue. Furthermore, this term is equal to half of the solar mass squared difference 
that is about 30 times as small as the atmospheric splitting (see table [T]). As a result, to 
recover a value for Am 2 tm close to the measured one, we need that the second term on the 
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right-hand side of eq. (39) is the dominant one. We can estimate the ratio between the 
two terms, calculating r, the ratio of the two mass squared differences: 

/ \-l _ ^ m atmNO,IO 

yr-Nojo) = ^—2 



Am sol 



2\ko\ 


|*i 


cos (9 ko - 9ki) A =F 




2 -2\k \ 


|*2 


cos(6» fco -9 k2 ))X 2 


4|* | 


|*i 


COS 1 


0k - 


9 kl 


)A 





(40) 



1 l'"J-l — I'-UM'"^! " \"«0 " re 2/^ 



2 



1*1 


2 


-2|feo 


|*2 


cos(# fco - fc2 ) 


4M 


fell 


COS 


(0fc o ~ ^l) 



The natural range of this quantity r _1 would be something like [0.3 — 0.7] (central value 
0.5 and corrections of order A). However, measurements give r = 0.0321qoo5 ~ A 2 , or in 
other words 1/r ~ 30 ~ A -2 . We conclude that 





|*i 


cos(6> fco - 9 kl ) 


|l*i 


2 . 


-2\ko\ 


*2 


cos (9 ko -9 k2 )\ 



X s . (41) 

The most natural explanation may be assuming that cos (9 ko — 9 kl ) is very small. In that 
case the absolute values of all parameters can still be of order one, which was part of 
the naturalness requirement of the model. Neutrinos present a quasi degenerate (QD) 
spectrum and both normal and inverse ordering are possible. The typical scale of neutrino 
masses vl is given in this case by 

«£ = -^u-^-o , Am ^ m , VI ~ (0-leV) 2 . (42) 

A 2 Wh] 2 - 2|ko||k 2 | cos (9 ko - 9 k2 ) 



Possible alternative solutions of eq. (41) that give a non-QD spectrum can be obtained 
only by admitting the parameters belong to a larger, but less natural, range, A 2 — A~ 2 . In 
this case when A; is of order A -2 and cos (9 ko — 9 kl ) ~ A while fci j2 still of order one a inverse 
hierarchical spectrum can be obtained. Another possibility to get an inverse hierarchical 
(IH) spectrum is having k and cos {9 ko — 9 kl ) of order one, k\ ~ A -1 and k 2 ~ A 2 . Finally 
the normal hierarchical (NH) spectrum can be obtained only in the case in which k 2 ~ A -2 , 
cos (9 ko — 9 kl ) and ki of order 1 and k ~ A 2 . 



3.2.3 Mixing matrices at Leading Order 

Apart from the charged fermion masses of the second and third families, mixing angles and 
masses at leading order do not fit the experimental data: the charged fermion first families 
are massless, the CKM matrix is the unity matrix and the PMNS matrix can account for 
two maximal rotations in the (12) and (23) sectors: 

/ l/y/2 -l/y/2 \ 
V = l, U = R 2B (~)R 12 (j)=\ 1/2 1/2 -l/y/2 . (43) 

4 4 V 1/2 1/2 +1/V2 / 

At this level of approximation, the PMNS matrix corresponds to the BM pattern discussed 
in the introduction. As already stated, this mixing scheme does not fit the data, due to 
the large value of the bimaximal solar angle which is out of about 5a. Only considering 
the NLO contributions, which we will study in the next section, the model agrees with the 
measurements. 
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3.3 Fermion mass matrices at higher orders 

The NLO contributions in the mass matrices originate from two sources: the first are the 
higher order terms in the superpotential, while the others come from the insertion of the 



NLO flavon VEVs in the operators in eqs. (22, 33). In section 3.4 we will show how the 



flavons develop VEVs and how they are corrected. Here we anticipate the results, reporting 
the flavon VEVs in a form which is useful for the discussion in this section: 



(<P) = I 1 j v v + ( | Sv v , {<ff) 
(X) = ( ] v x + [ 1 | 6 2 v x , (a) = v a . 













v + 1 











\0J 





(44) 



Some comments are noteworthy: the subleading corrections are suppressed with respect 
the LO terms as 5v/v ~ A and S 2 v/v ~ A 2 for each flavon; NLO corrections to the second 
and third entries of (ip) and (<£>') are present, but they present the same structure of the 
LO terms and can be re-absorbed; similarly, the NLO corrections to the third entry of (%) 
and to (er) are present, but they can be re-absorbed into the LO terms; the other entries of 
(x) do not receive any corrections at NLO. If we consider the NNLO approximation level, 
i.e. corrections of relative order A 2 with respect the LO terms, we see that the second and 
the third entries of (</?) ((<p')) are not (anti-) aligned anymore and that the second entry of 
(x) is filled in. It is interesting that the first entry of (%) is still vanishing at this level. We 
will see the relevance of this structure in a while. 

3.3.1 Dirac mass terms 

The Dirac matter superpotential at NLO is given by 

3 1 
E2/6,(^FlF 3 c (0 + , )^ (5) + 



A 2 

i=l 



3 8 

+ E ma ^FlF^pxP + Y J y m ^,F L F c 2 e\<p + <fix t 



i=l ' i=l 

6 7 



, (?) + 



(45) 



i=l ' i=l 



i=l 
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where, as for eq. (p2|), we adopt a compact notation: 



X 



(5 



X 



(6 



(7 



X 



(8 



(9 



{y? 3 x, vV'x, </V 2 x, y? /3 x, v 3 ^, vV^> </W v'M 
{y? 2 x, w'x, ^' 2 x, yV, v'M > 



(46) 



(10) 



<^'V x 5 , xV x 3 ^ 2 } • 

Note that not all of these terms are non-vanishing when the flavons develop VEV: in 
particular Xf 1 ' for any % do not give a contribution when the LO VEVs are considered: only 
when the corrections to the VEVs are introduced, they contribute to the mass matrices. 
For this reason also the terms with X^ must be taken into account, even if they are 
suppressed by an additional A. The other terms which are vanishing at this order of 
approximation are xj; for i = 4,5,6 and X^ for i — 6, . . . , 13. When flavons and Higgs 
fields develop VEVs, we get the following Dirac mass matrices: 



M. 



NLO 



-3 



Af 



NLO 



Al 



NLO 





m o 

V9 o 



m o 

m o o 



m o o 

yio o o 



v d p X 5 - 3 



v d p X 5 + 



v u A 7 + 






m A 











\ 





2/2 


















~2/2 
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^A 






' 





2/eA \ 





fa. 


I) 


v d p X 3 + | 













-fa 











Vi J 





m A 


N 









2/eA ' 





2/3 





J ^A 5 + 


(i 





fa 





-fa 


o , 









fa , 



<A ,(47) 



v}\ , (48) 



^A,(49) 
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and m^) LO = M^ LO , where we used the definitions 



2/i 

2/2 

m 

2/4 

m 

2/6 
2/7 
2/8 
2/9 

Vio 



2/i +2/6,(i)A , 
2/2 + 2/7,(i) A , 

1 



9,(i) 



A 



fay fa^ 



^2 [2/ll,(0] + ^3 



2/io,(i) : 



fay fay 



1 fa, 



A 2 



2/5; 



5V 



(2/6,(2) +2/6,(3)) + v2/l 



P 

Uy 

fay 



(50) 



(2/7,(2) +2/7,(3)) + v2/2 

A ty 

Z^i=l 2/8,(0 + T zU=2 2/3,(0 



ty 



^5 [2/ 



9,(0 



"2/4 



fay _^ fay 



•^6 [2/U,(0] 



^7 



2/io,(i); 



fay fay' 

f y ' v 



In the previous definitions we can see that each y is the sum of two pieces: the first refers 



to the terms in eq. (45) when the LO flavon VEVs are considered; the second comes from 



the terms in eq. (22) where the NLO flavon VEVs are introduced. The only exception is 
J-" 3 which refers to the term proportional to X-- in eq. (45) and that give contribution 



only when the NLO flavon VEVs are considered. These two parts are of the same order of 
magnitude, since 5v/v ~ A and 5 2 v/v ~ A 2 . Note that Ti are distinct linear combinations 



of the arguments in the squared brackets. The expressions in eqs. (47, 48, 49) are valid 



at NLO level and note that the ( ant i-) alignment between the second and third entries of 
the (second) third families are still preserved. When considering higher order terms, this 
feature is lost and the (1, 1) entry of each mass matrix is filled in. 



The values for the charged fermion masses given in eq. (28) are modified only by sub- 
stituting the coefficients yi with their tilde-versions: 



2/1^2/1, 



2/2-^2/2, 



2/3-^2/3- 



(51) 



At this approximation level, the first family masses are not yet well-described, because 
they are too small. 
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3.3.2 Majorana mass terms 

Moving to the Majorana part of the matter superpotential at NLO, we get 



A 3 

i=i 



+ ^F^A R + 



+ ^ W 2 A W + E ^ W 2 A fl x£? + E fW 2 Arf?+ 

i=2 i=4 

+ E ^ W 2 A*Af> + f ^ w 2 A«x£t>+ 

i=l i=5 
Z 5 mr. 77C/)4 , 



+ E Z ~^F?F^X^ , 



i=i 



(52) 



where as usual we used the compact notation with 

Af x) = {x 3 ,x 2 ^ 3 }, 
x f 2) = Wx, , 

xf 3) W'X 2 , ff V} , (53) 

^ (14) = W 3 x, y?Vx, vV 2 x, y? ,3 x, <p 3 °, y?V^, vW v>'M , 

A few comments are in place. Note that all the terms proportional to k 3 can be reabsorbed 
by a redefinition of &o,i,2 anc ^ that the only new structure which corrects Ml comes from 
the term FlFlAlX when we consider the subleading corrections of (%): as a result the 
entries (1,3) and (3,1) of Ml are filled in by terms proportional to A 3 . Regarding the 
contributions to the Majorana mass matrix for the RH neutrinos, it is important to say 
that the terms proportional to z 3 ^ with i = 1, . . . , 3 and to z^ with i = 1, ... ,4 are 
vanishing, due to the particular flavon VEV alignment of the model. As a consequence all 
the NLO contributions to Mr are of the order of A 6 , apart that one to the entry (2,2) 
which is of the order of A 4 . Finally, we note that each entry of Mr is independent from 
all the others, being Ff singlets of the flavour symmetry, and therefore all the Zi are free 
parameters with modulus of order 1. We listed only the dominant contributions, but the 
higher order terms would correspond to subleading corrections, which we can safely neglect 
in the following. 
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As a result of this analysis the Majorana masses for LH and RH neutrinos are given by 



k' k[X k' 3 X 3 \ ( z 6 X 6 z 5 X 6 z 4 A 6 

M N L LO = ( k[X k' )v L , M% LO = z 5 X 6 z 2 A 4 z 3 A 6 I v R , (54) 
^A 3 k' + k' 2 X 2 J \ z 4 A 6 z 3 X 6 Zl 

where with the notation k[ we account for all the redefinitions done on the parameters. As 
already stated when discussing the LO mass matrices, the contributions to the effective 
light neutrino mass matrix come from the type-I and type-II See-Saw mechanisms. The 
resulting NLO m[/ type is given by 

m (*yp^ =m N D LO {M^ LO )-\m N D LO ) T 

vl>? MeA yiy^ \ ( v ^yx 2 (55) 

yiyeX yl yi — , 

MeA yl yl j ZiVr 

and it is diagonalized by a maximal rotation in the (23) sector and not in the (12) sector 
as demanded by eq. As a result, we need that the type-I See-Saw contribution is at 



least 0(A ) the type-II See-Saw one, which we can express using eq. (51) as 

^- < X%v L . (56) 

ZlVR 

We remember here that vl,vr are the VEVs of Al and A# respectively. In particular vl 
is the VEV developed by the SM (1, 3, 1) triplet component of and it is induced once 
the EW symmetry is broken. As we will see in detail in the next sections the physical SM 
triplet T ~ (1,3, 1) arises by the mixing between the SM (1,3, 1) components of A L and 
of two additional fields, S and £' transforming under the PS gauge symmetry as (1,3,3). 
The VEV (T) of the SU(2)l triplet T is related to its mass Mt through the following 
expression 

where are numerical coefficients arising by the details of the scalar potential and v\ = 
(h u ) , v 2 = (h' u ) are the VEVs of the two up-type light Higgs doublets needed for the 
realization of our model (see appendix |B| for details). Since vl is the projection of (T) 
along A L neglecting fine tuned cases it holds that 

v L ~ (T) . (58) 



From eq. (56) we need therefore that Mt and vr satisfy 



M T < X [ — 2 3 - ) v R . (59) 



At the same time neutrino mass data imply that 



k' 0aij ^<O(l)eV. (60) 
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Combining the constraints of eqs. (59)-(60) we see that in the most natural scenario, 
assuming Oy = 1, My and vr satisfy 



v R >30M T , 

(61) 

10 12 GeV < M T < 10 13 GeV . 

Nevertheless if we allow the numerical factors laying in the range 0.1 — 10, then M?, 
and consequently vr, can be reduced even of two and one orders of magnitude respectively 

10 10 GeV < M T < 10 12 GeV . (62) 

In the following discussion of lepton mixing and in the phenomenological analysis, we 
will assume that indeed type-II See-Saw is dominating and we will neglect the type-I 
contributions. In the section devoted to the study of the scalar potential we will justify 
and find out the region of the parameters space where indeed type-II See-Saw dominates 
over type-I. 



3.3.3 Mixing angles at the Next-to-Leading Order 



Looking at eqs. (47)-(49)-(54) we see that the fermion mass matrices are of the required 
form as in eqs. ( 10 )— ( 13 ). The resulting mixing matrices are modified with respect to 
the LO approximation and interesting new features follow. On the quark sector, the 
CKM matrix receives deviations from the unity and at NLO the angle 9\ 2 is not vanishing 
anymore: 

Q q * 13/73/3 -2/22/8 1 /ftQA 

V2 1 2/2 1 m\ 

Looking at this result, the meaning of the parameter A is clear: it is defined as the ratio of 
the flavon VEVs over the cut-off of the theory, but it also determines the order of magnitude 
of the Cabibbo angle. This justify our initial assumption of A = 0.2. We note that this 



result is in concordance with the discussion below equation (13). If the second columns of 



the up- and down-quark matrices are not proportional to each other, we can generate a 
non- vanishing Cabibbo angle 9f 2 , while the two other angles in the CKM matrix are still 
vanishing. 

In the lepton sector, the PMNS matrix corresponds to the BM scheme which receives 
the corrections as illustrated in eq. @ and as a result the solar and reactor angles are 
given by 



13 



12 



A 
2 

7T 

4 



2/6 

2/1 



2/7 
2/2 



2/i 2/2 



(64) 
(65) 



As it is easy to see, the reactor angle and the deviation from the maximal value of the solar 
angle are of order A and therefore the model is now in agreement with the experimental 
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dataj^J This corresponds to a weaker version of eq. ([!]). At this approximation level, 
the atmospheric angle remains maximal. As already stated, the forthcoming neutrino 
appearance experiments will provide a good test for the model, indeed they will be able to 
verify if 9[ 3 is close to its present upper bound. 



3.3.4 Higher order effects 



At the next-to-next-to-leading order (NNLO) and even higher orders, many new terms 
appear in the superpotential. However, only few of them lead to new terms in the mass 
matrices, while the rest can be absorbed in redefinitions of the parameters as in eq. (50). 



For this reason we do not report the full list of NNLO contributions, but we just comment 
on the physical consequences. Three effects are worth mentioning: 

- as expected, the masses of the first families are strongly suppressed. We find that 
the down-quarks and the electron masses are suppressed by a factor of A 6 and the 
up-quark mass by a factor of A 8 . This leads to the following mass hierarchies: 







m e 















m,. 



m t 



A 7 



(66) 



As it is easy to see, the electron and the up-quark masses perfectly fit the experimental 
values and there is only a small tension in the down-quark sector, where we expected 
a smaller suppression. However, the large number of parameters, in particular from 
0, 0' as well as p couplings, which enter in the definition of the first families allows 
to correctly fit the down mass, without affecting other observables. 

the (anti-) alignment of the (23) and (33) elements of the Dirac mass matrices in eqs. 
(47-49) gets broken as the new terms appear. The new elements are A 2 suppressed 



with respect to the older terms. As a result, the matrix that diagonalizes MjM/, 
with % — e,u, d, has no longer an exact maximal mixing in the (23) sector. In the 
lepton sector, this translates to a A 2 deviation from maximality in the atmospheric 
angle of the PMNS matrix. In the quark sector, the angle #23 becomes of order A 2 . 
It is interesting to note that 9\ z remains vanishing at this order. It only appears 
when even stronger suppressed terms are taken into account and is of order A 3 , in 



accordance with the Wolfenstein parametrization 31 



the third columns of the mass matrices in eqs. (47)-(48) are proportional to v^, 
while the second column of M^ LO is proportional to fjand that one of M^ LO to 
— 3Vp. Therefore, also at NLO, eqs. (p0| are fulfilled. At the NNLO level, terms 
proportional to v d p appear in the thirdcolumns of charged lepton and down-type 
quark matrices and terms proportional to v d i in the second columns. The new terms 
are A 2 m 5% suppressed with respect to the old entries. We thus expect deviations 



from the relations 



m. 



\m b 



and 



m, 



3|mJ at the 5% level. 



5 Deriving eq. (65 1, we already neglect the corrections which increase the value of the solar angle, instead 
of decreasing it 



22 



3.4 Flavon scalar potential 

In this section we comment on the vacuum alignment mechanism which explains the flavon 



VEVs as in eqs. ([17)— ( 18 )— ( 19 ) . It turns out that our desired alignment is exactly the 
one presented in |23| , once we transform all the fields in our basis. Notice indeed that it 
is possible to identify each flavon of table [3] with the flavons in 23 , by simply comparing 
the transformation properties under the full flavour group: 



9 



Xi 



X 



a 



Using the following unitary matrix to move from the basis in |23| to our basis, 

10 

-i/y/2 i/y/2 
1/V2 l/y/2 



(67) 



(6* 



we find (up to irrelevant phases) the following flavon VEV alignment]^] We will comment 
in a while about the presence of equivalent solutions. 
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(69) 



which correspond to eqs. ( 17 )— ( 18 ) . 

In [23], a set of driving superfields has been introduced to guarantee the correct flavon 
vacuum alignment: these new fields are gauge singlets and transform only under the flavour 
group, but differently from the flavons, they do not develop vacuum expectation values^} 
Under the continuous i?-symmetry all the driving fields transform as U(1)r = 2 and 
therefore, constructing the superpotential, they appear only linearly. For the same purpose, 

In table |5l we 



we introduce in our model a set of driving fields which recall those in 23 
show the driving fields and their transformation properties under S4 x Z4. 

It is easy to determine the correspondence between our set of driving fields and those 



of 23 



D, 



<Pr 



Xb. 



ob 



(70) 



We now construct the driving superpotential Wd, which contains only flavons and driving 
fields and in particular neither matter fields nor Higgses, and look for the conditions that 
minimize the scalar potential, 



+ m 2 |$ i | 2 + 



(71) 



6 The VEVs of the fields ip and ip 1 is recovered by applying the unitary matrix in eq. ( 68 ) to an equivalent 
configuration of eq. (51) in 23 , resulting from the application of the element (TS) 2 to eq. (51) and of the 



element T to the eq. (18) of [23] 

7 This strictly holds only in the exact supersymmetric phase, while in the broken phase they develop a 
VEV proportional to the common soft breaking scale [32], usually denoted as msuSY- This could have a 
relevant impact when discussing on flavour violating processes, as studied in a series of paper [33| dealing 
with a context similar to our model. 
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Table 5: The driving field content and their transformation properties under S4 x Z4. They are all singlets 
under the gauge group and the FN symmetry, while they transform as U(1)r — 2 under the continuous 
R-symmetry. 



where $j denote collectively all the scalar fields of the theory, m? are soft masses and dots 
stand for D-terms for the fields charged under gauge group and possible additional soft 
breaking terms. Since ntf are expected to be much smaller that the mass scales involved 
in Wd, it is reasonable to minimize V in the supersymmetric limit and to account for soft 
breaking effects subsequently. 

we 



Since our flavon and driving field content exactly corresponds to the one in 23 



already know that the VEV alignment in eqs. ( 17)— ( 18 ) represents an isolate minimum of 
the scalar potential. We only need to identify the relations which link the VEVs v v , ty, v x 
and v a among each other in our model. To this purpose we write the driving superpotential, 
which reads as 



w d = fiD R (f<f + f 2 D R (p'(p' + f 3 D R <f<f' + U<p R (f(p'+ 
+ M^xrX + fhXRXO + feXRXX+ 
+ M\ A 2 a R + M 3 Aa R a + f 7 a R aa + f^RXX , 



(72) 



where the first line deals with only the fields ip and </?' and the other two with x and a. As 
also explained in [23], this leads to the alignment in eqs. (17, 18) where the VEVs satisfy 
to 



fiv% + hify + v^/s^VV = 



Mi 



Vlh 



(73) 



The solution in eqs. ( 17. 18 . 73 ) is not unique, but it is possible to introduce a set of soft 



supersymmetric breaking parameters, which selects this solution as the lowest minimum 
of the scalar potential. 

It is interesting to note the presence of an other source of uncertainty in our solution, 
which minimizes V. Given the symmetry of Wd and the field configurations of eqs. (17, 
18, 73), by acting on them with elements of the flavour symmetry group S4 x Z 4 , we can 
generate other minima of the scalar potential. These alternative solutions however are 
physically equivalent to those of the original set and it is not restrictive to analyze the 



model by choosing as local minimum that one in eqs. (17, 18 73). 



For the FN field 9 to acquire a VEV, we assume that the symmetry U(1)fn is gauged 
such that 9 gets its VEV through a D-term. The corresponding potential is of the form: 

1 



V, 



D,FN 



-(M FI - g FN \9\ 2 + 



(74) 
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where gpN is the gauge coupling constant of U(1)fn and M FI denotes the contribution of 
the Fayet-Iliopoulos (FI) term. Dots in eq. Q74j) represent e.g. terms involving the fields 
F{ and F% which are charged under U(1)fn- These terms are however not relevant to 
calculate the VEV of the FN field and we omit them in the present discussion. Vd,fat leads 
in the supersymmetric limit to 

M 2 = K*)| a = ^. (75) 
9fn 



It is relevant to underline that the VEVs in eqs. (73 75), depend on Mass param- 
eters: all these mass scales naturally have the same order of magnitude and as a result 
VEVs/ Af ~ A. The only exceptions are the VEVs of (p and (p f , which depend on a flat 
direction. In the model, we simply assume that their VEVs have values of the same order 
of all the other flavon VEVs. 

3.4.1 Higher order contributions 

In this section we briefly comment on the corrections which enter in the flavon VEVs, 
once the higher order contributions are taken into account. We leave all the details to the 
appendix [Cj 

In the superpotential Wd, the flavons which contribute to the Dirac mass terms, ip 
and tp', and those which contribute to the Majorana mass terms, x an d c, at LO belong 
to two separated sectors, indeed any mixing term is prevented due to the Z 4 symmetry. 
This situation is not preserved at NLO, since the fields \ and a are neutral under the Z4 
symmetry and therefore we can add each of them to all the terms in Wd- This leads to 



modifications to the LO VEV alignment of eq. ( 17) and it turns out that the first entries of 
(<p) and (<p') are filled in, while the second and third entries are corrected by terms which 
can be however absorbed into the LO ones, without spoiling the alignment. Also the VEVs 



in eq. (18) receive some corrections: the first and second entries of (x) still vanish and 



the NLO contributions to the third entry can again be absorbed into the LO term. This 



discussion justifies the results showed in eq. (44). 



Corrections from the NNLO contributions are without particular alignments: in par- 
ticular, the second and third entries of (<p) and (tp') are no longer related, and also the 
second entry of (x) gets a non-zero value. It is interesting to note the the first entry of (x) 
remains zero. 



4 Higgs scalar potential 

In this section we present the study of the Higgs potential in our model. It is an interesting 
example of how the introduction of flavour symmetries and the assumptions done to get 
the correct mass matrices have non-negligible consequences on the Higgs sector. As a 
result, the study of the Higgs scalar potential and of the gauge and the Yukawa coupling 
runnings in a general non-flavour PS context does not strictly hold. Notice that even if the 
following analysis refers to our particular choice of fields and symmetries, our conclusions 
can be taken as a general hint for a very large class of models that combine a discrete 
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flavour symmetry with a grand unified scenario: indeed our model building strategy shares 
common features with other constructions. In particular the Higgs fields usually transform 
under the flavour symmetry Gf and this has direct consequences on the implementation of 
the grand unified symmetry breaking. Moreover type-II See-Saw dominance and particular 
patterns of vanishing projections of the heavy Higgs fields on the light Higgs doublets are 
frequently required to get the correct fermion mass matrices. 

In table [6] we list all the Higgs fields which are necessary to reproduce the correct mass 
matrices and to implement the desired PS symmetry breaking pattern: in the first part 
of the table we report the Higgses already introduced in table [4], while the second part 
contains all the additional fields. 



Higgses 


<f>,# 


P 


A L 




PS 

z 4 


(1,2,2) 
1 


(15,2,2) 
-1 


(TO, 3,1) 
1 


(10,1,3) 
-1 



Higgses 


A L 




A 


B 






i 


PS 

z 4 


(10,3,1) 
1 


(10,1,3) 
-1 


(15,1,1) 
1 


(15,1,1) 
-1 


(1,3,3) 
1 


(1,3,3) 
-1 


(i,i,i) 
-i 



Table 6: All the Higgs fields of the model and their transformation properties under the gauge group and 
under Z±. Notice that they are invariant under the other factors of the full flavour symmetry group Gf. 



Some of these Higgs fields are already present in the minimal version of PS models [34] : 
typically a (15,1,1) multiplet - as the fields A and B in table [6]- is used to break SU(4)c 
down to SU(3) C x U(1) B -l and to induce the VEVs of the couple of (10, 1, 3) © (TO, 1, 3) - 
corresponding to the fields A R (BA R in table [6] The latter VEVs breaks SU (2) R x U (1) b-l 
into the SM hypercharge U(l)y concluding the symmetry breaking chain from the PS gauge 
group to the SM SU(3) C x SU(2) L x U(l) Y - The field that triggers the EW symmetry 
breaking is usually a bidoublet(l,2,2) - as or </>' in table |6j The fields (10, 3, 1) © (10, 3, 1) 
- the fields A^ © A^ in table |6| - do not develop VEVs at tree level in the usual minimal 
PS model, but only when next to leading order terms are taken into account; these are 
typically suppressed by the Planck scale, as already stated in [34]. For this reason in the 
minimal PS the type-II See-Saw contributions to the effective neutrino masses are almost 
negligible. 

We identify three main reasons for which the existent studies of the symmetry breaking 



patterns in the PS context 24,34 have to be modified and this automatically justifies the 



presence of the new fields in table [6] 

- the assumption v u p = necessary to distinguish the up-quark sector from the others 
and to recover the up-quark mass hierarchies can be realized only if we include two 
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identical copies of bidoublet (1,2,2), an d 0' (see details in appendix [B]) , and we 
then impose that four SU(2)l doublets (2 up-type and 2 down-type) remain light; 

- since the fields p, Ar and Ar transform non-trivially under the flavour symmetry Z±, 
it is necessary to introduce two copies of (15,1,1) multiplets, A and B, with opposite 
Z4 charges, 1 and -1 respectively: A is responsible of inducing the breaking of SU(2)r 
through its coupling with Ar and A^; B allows the coupling of the bidoublets 0, 0' 
with the bidoublet p. In this way all of these three fields have a non-vanishing 
projection on the light Higgs SU{2) l doublets; 

- the component of A L which corresponds to the usual SM triplet (1,3,1) can develop 
a VEV once the EW symmetry is broken only in the presence of a trilinear coupling 
with the SU(2)l Higgs doublets. This coupling cannot be originated with only the 
Higgs fields A^ © Al and the field content given in table [6] For this reason we need 
an additional field which mediate this coupling and the simplest choice would be a 
bitriplet E ~ (1, 3, 3) that can couple with the fields 0, 0' or p and at the same time 
can mixes with A^, when Ar develops VEV at the SU(2)r breaking scale. However, 
once more, the presence of the Z4 symmetry obliges the introduction of two distinct 
(1,3,3) Higgs fields, E and E', with opposite Z 4 charges, 1 and —1 respectively. In 
this way E can couple to the bidoublets 0, 0' or p, while E' can mix with A^. Finally, 
we need a new ingredient that mixes E with E': a PS singlet £ charged —1 under Z4 
can do the job. 

The scalar part of the superpotential is then given bj 
W = 00 + X -My 0'0' + M H , 00' + \m p pp + M Al A l A l + 

+ M Ar A r A r + l -M A AA + X -M B BB + ^M S EE + ^M^E'E' + ^M&+ 

+ A C EE'£ + XzabABS, + X^Bp + \^ p <p'Bp + \\ A AAA+ 



3 

+ l -\ B BBA + \ L A L A~ L A + \rA r A~rA+ 

+ ^A < j !)E 00E + ^A0/ S 0'0'E + A^ E 0'0E + ^A pE ppE+ 

+ Aas'A^AkE' + Aas'A^E' + ^A E EEE + ^E'E'E . 



(76) 



The vacuum configuration at the GUT breaking scale is given b)|j 

(A R ) = (A R ) = Mr , (A)=M Cl , (B) = M G2 , (0 = V ( . (77) 



5 Since all the Higgs fields are neutral under the continuous U(1)r, the scalar superpotential explicitly 



breaks it, while preserving the usual R-parity. The terms in eq. (76) could be generated from a U(1)r- 
conserving superpotential in which the breaking is mediated by additional fields, which are U(1)r = 2 
and develop non-vanishing VEVs. For instance the mass term (jxf> could originate from a trilinear term 
X (jxfi, when (X) — M^. Similarly the trilinear coupling A^SS'^ could originate by the non-renormalizable 
term X££'£/A', when (X)/A' = A^ and A' is the energy scale of the dynamics of the field X. In our 



model we simply assume the existence of the terms in eq. ( 76 1 in the superpotential and allow for an 
explicit breaking of the U(1)r symmetry in this sector. 

9 We redefine (A^) = vr as (A^) = Mr in order to adapt to the usual notation. 
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The VEVs (A) and (B) break the SU{A) C to SU(3) C x £/(1)b-l, while (A R ) and (A R ) 
break SU (4)c x £{7 (2)# into S*?/ (3)c x Z7(l)y. Therefore given our field content the colour 
breaking scale, Mc = Max(Mc 1 , Mc 2 ), is never smaller than the SU(2) R breaking scale, 
Mr. We may expect that M Cl ~ M C2 = M c but in principle they could be different. 
Finally is expected to be close to the flavour breaking scale, due to the £ gauge singlet 
nature. The PS breaking pattern so far sketched is therefore summarized by 

SU(4) C x SU(2) L x SU{2) R 
I 

SU(3) C x U(1) B -l x SU(2) L x SU(2) R (78) 
I 

5'[/(3) c x SU(2) L x C/(l) y . 



The F-derivative system obtained by the superpotential in eq. (76) reads as 



m Ar + ^rM Ci = o , 

M B M C2 V2X B M Cl M C2 + \ AB ^M Cl V^ = , 

M A M Cl + ^|AbM£ 2 - ^R m r + V2X A M^ + \abzM C2 V(: = 
MM + X ABi M Cl M C2 = . 



(79) 



By solving the previous equations, we can express the mass parameters that enter in 
the superpotential in term of the adimensional parameters A« and the physical breaking 
scales. All the details regarding the mass spectrum are reported in the appendix |Bj but 
some comments are in place. As in the minimal supersymmetric PS |34| when the singlet 
component of A develops a VEV, there is an accidental SU(3) symmetry involving A R and 
Ar. When the singlet components of these fields acquire a VEV the accidental symmetry 
is broken to SU{2) giving rise to 5 Goldstone Bosons (GBs). At the same time SU{2) R x 
U(1)b-l is broken down to U(1)y, eating 3 of the 5 GBs. Therefore 2 of them, namely 
5 ++ and 5 are left massless, down to the SUSY soft breaking scale ~ 1 TeV. This is 



a well known prediction of SUSY PS theories, which can be tested at LHC 35 . On the 
other hand, contrary to the minimal case described in [34] , due to the mixing between the 
Higgs fields A and B, no colour octet is lighter than M R . 

In order to assure type-II dominance and to get the correct PS symmetry breaking 
pattern, it is necessary that Mj- < M R < Mc < V%. We will see in a while the constraints 
on the parameters which can be derived from this special mass ordering. For the moment, 
we just assume this scheme and we describe the mass spectrum in terms of the SM gauge 
group. Starting from the heaviest energy scale 

- at V f 



2 heavy SU(2)l triplets given in sec. B.3.2 
between and Mc 
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B.1.2 



B.1.3 



B.l.l except one, called £o, 



B.2.1 



all the SM singlets given in sec. 
the colour triplets given in sec. 
the colour octets given in sec. 
the heavy doublets given in sec. 
the two heavy couples of 577(2) l triplets given in sec. 

at M c 

the color scalars originating by A R © A R , p and A L 



B.3.1 



the SU(2)l singlets given from sec. B.1.6 



to sec. 



B.1.10 



the SU(2) L doublets given from sec. B.2.2 



to sec. 



B.2.4 



the SU(2) L triplets given in sec. B.3.3 and sec. B.3.4 



at M R 

the singlet £o! 
at M T 



the light couple of SU(2)l triplets given in sec. B.3.1 
at Msusy 



the scalar singlets 5 ++ and S given i n sec. 
the SU{2) L light doublets given in sec. B.2.1 



B.1.5 



With the Higgs field content given in tables [4] and [6] and the scalar spectrum so far sketched 
we can plug the gauge coupling runnings and see if the conditions given in eq. (61 ) or those 



in eq. (62) can be satisfied. Furthermore, in the study of the gauge coupling runnings (see 



appendix D for details) from the Mqut to the EW scale we have to impose the following 



constraints: 

- recover the EW values for 03,02 and a±, related to the gauge couplings of the SM 
gauge group ; 

- impose that, at Mr < Mc, U(1)y is originated by the SU{2)r x U(1)b-l breaking; 

- impose chr^l = etc at Mc ; 

- define the GUT scale as the scale at which the largest = 1. In this way we are 
sure to be in a perturbative regime up to the GUT scale and thus we are allowed to 
adopt the one-loop renormalization group equations (RGEs). 

Even using all these constraints we are left with two more freedoms, the value of the S77(4)c 
and SU{2)r breaking scales, i.e. Mc and Mr respectively. 

We adopt two distinct approaches, that we will indicate as the more constraining and 
the less constraining ones. In the first case we define Mc the scale at which the largest a.i 
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is equal or smaller than 1/4-7T. In this way all the gauge coupling at Mc are smaller than 1. 
In the second case we allow the largest oti to correspond to a gauge coupling in the range 



1 < <7i < 3. Then Mr should satisfies eq. (61) or eq. (62), but its exact value is not fixed 
yet. 

A few general comments are needed. The non-minimal PS field content affects the 
gauge coupling runnings in a non-negligible way. In particular the presence of the charged 
singlets 5 ++ and 5 down to Msusy deeply modifies the U (l)y and SU{2)r gauge coupling 
evolution. It turns out that the largest above the Mr scale is always olr. Therefore the 
two approaches we described can be formulated as follows: 

- More constraining approach -<=>- olr < 1/47T, at Mc 

- Less constraining approach -<=>- 1 < qr < 3, at Mc ■ 

More constraining approach 

In this case there are no solutions neither for the ranges of values of Mt and Mr given 



in eq. ( |6T| ) nor for those in eq. (62). Indeed we find that Mr < 10 12 GeV, as can be 
seen in fig. [TJ In other words if we adopt this constraining approach to fix the value of 
Mc, the type-I and type-II See-Saw scales require Yukawa parameters which are at least 
2 order of magnitude far from their natural values to reproduce the correct neutrino mass 
scale. We cannot be satisfied with such a solution because in this case the type-II See-Saw 
dominance is obtained by increasing the Yukawa couplings in the right neutrino sector and 
at the same time reducing the coupling of the left-handed neutrinos with the scalar triplet. 
Even if this may be considered a solution, our challenge was to provide a justification 
of type-II See-Saw dominance through the analysis of the Higgs scalar potential and not 
by tuning the Yukawa parameters. Moreover we introduced a FN Abelian symmetry to 
explain the small (< 10 -2 ) Yukawa parameters necessary in the charged fermion sector to 
reproduce the correct mass hierarchies. The presence of Yukawa parameters of this order 
in the purely left-handed neutrino sector makes the introduction of the FN symmetry 
questionable. Notice that in the discussion of the neutrino spectrum we have enlarged the 
parameter range up to values not smaller than A 2 . 

Less constraining approach 

In this case there are solutions only for the second range of values of M? and Mr given 



in eq. (62). As can be seen in fig. |2j Mr can now reach the value of 10 13 GeV. However 
the three scales Mr, Mc and Mcut are compressed in a narrow region around 10 13 GeV 
and therefore our model is described by an extended MSSM model almost up to the GUT 



scak-^ Nevertheless, its PS origin is reflected in the non-trivial relations between the 
Yukawa couplings. In conclusion, by admitting the Yukawa parameters span in a range 
0.1-10 and switching Mr, Mc and Mcut very close to each other, we have found a narrow 



10 Above this energy scale an other gauge structure could be active. We do not take in consideration an 
high-energy completion of the model, but it is reasonable that larger gauge groups or particular construc- 
tions could be present at these energies: for example an SO{lQ) inspired approach in which fermions do 
not belong to a unique representation. 
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Figure 1: The running of the gauge coupling constants in the more constraining approach. Mt = 10 11 
GeV, M R = 10 12 GeV, M c = 1.4 x 10 12 GeV (where a R = 1/Att) and M gut = 4.0 x 10 12 GeV (where 
ctji = 1). In the dotted figure, we show a detail of the full plot, restricting the energy scale inside the range 
10 io^ 10 i3 Q e y. 



region of the parameter space where our model could still give a realistic description of 
fermion masses and mixings and in which type-II See-Saw dominance is not imposed by 
hand. 

To finally consider our model viable, we should study the stability of the flavour struc- 
ture of the mass matrices under the RGEs from the GUT scale down to the EW one. The 
study of the full set of the RGEs of the model presented is beyond the purpose of this 
paper. For this reason we will neither run the parameters of the scalar superpotential nor 
include and run the parameters of the soft SUSY breaking potential. Under these approx- 
imations the EW vacuum expectation values do not change from the GUT scale down to 
the EW one. However this does not affect our conclusions for what concerns the stability of 
the mass matrix structures since the EW VEV shifts due to the running factorize out and 
leave the Yukawa flavour structure unchanged. The study of the Yukawa matrix running 
will be done in the next section. 
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10 2 10 4 10 6 10 8 10 10 10 12 10 14 

fi (GeV) 

Figure 2: The running of the gauge coupling constants in the less constraining approach. Mt = 10 12 GeV, 
M R = 10 13 GeV, M c = 1.8 x 10 13 GeV (where a R = l/3) and M gut = 2.2 x 10 13 GeV (where a R = l). In 
the dotted figure, we show a detail of the full plot, restricting the energy scale inside the range 10 11 -f- 10 14 
GeV. 

5 Yukawa coupling running 

In the previous section we analyzed the constraints on the scalar Higgs sector coming 
from the requirement of type-II See-Saw dominance and from the presence of the flavour 
symmetry under which the Higgs fields non-trivially transform. We found that the model 
is viable only in a small region of the parameter space for which Mr, Mc and Mgut ~ 10 13 
GeV are very close to each other. At the same time Mt lies at only one order of magnitude 
below Mr- For these reasons to study the stability of the flavour structure of the fermion 
mass matrices at low scale we can consider only the running from Mt onwards, thus 
neglecting the running from higher energies. Furthermore we also neglect the running from 
Msusy to the EW scale, which would introduce only minor corrections. We work under the 
assumption that type-II See-Saw is dominating over type-I and moreover that the effects 
from the type-I terms under the RGEs are negligible. Therefore in the studying of the 
Yukawa coupling running we do not take into account the Weinberg operator originating 
by integrating out the right-handed neutrinos. The error we introduce in this way is less 
than A 2 and we will see in a while that these contributions do not modify our results. 
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Furthermore, we study the stability under the Renormalization Group (RG) running, in 
the approximation corresponding to the NLO, i.e. considering the mass matrices introduced 
in eqs. p7b-(|49b. 



5.1 Yukawa matrices at Mt 

Since we start the renormalization group running at Mt, at this scale we integrate out the 
SU (2) l scalar triplet T obtaining an effective Weinberg operator responsible of the type-II 
See-Saw contribution. We recall here the origin of this effective operator. The Majorana 
parts of the matter superpotential given in eqs. (33), (52) contain terms with the coupling 

F L F L A L , (80) 



while the scalar part of the superpotential in eq. ( 76 ) contains the terms 

2 _ _ _ i i ( 81 ) 

+A AE ,A L A R E' + A A e'A r A l £' + -A S EEE + -A^E'E'E , 

that ensures the mixing between the (1,3,1) ((1,3,-1)) components of A^ (A^), E and E', 
whose lighter combination is identified with T (T), and provides the coupling of T (T) 
with the light doublets h d and h' d (h u and h' u ). The effective Weinberg operator at Mt is 
given by 

T T h h 

^Y Lrs M ; (82) 

where is the SU(2)i lepton doublets, h Ul = h u , h u , 2 = h' u and a^- are coefficients arising 
by the scalar potential and Y L is given by 

/ k' k[X \ 
Y L = k[X k' . (83) 

\ k' + k' 2 X 2 J 

Notice that we are neglecting the A 3 terms in Yl, because they would be irrelevant for the 
following analysis. 

For what concerns the charged fermion Yukawa, at Mt the Dirac part of the superpo- 
tential is written as 

Y U Q U c h u + Y' U Q U c h' u + Y d Q D c h d + Y' A Q D c h' d + Y e L E c h d + Y' e L E c h' d , (84) 

with 
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/ y 8 X 5 y 6 X \ 

1 \yi^ -y 3 A 4 yi J 



( y 7 X 

Y d = U 13 y A X 2 m 

\ Z/gA 2 -fa 

I y 7 X 

Yd = U23 y 4 A 2 y 2 

\ ygX 2 -y 2 



\ / y 6 X \ 

A 2 + yi , 
/ V m J 

\ / y 6 A \ 

A 2 + P , (85) 

/ V Vi I 



/ y 7 A \ / y 6 X \ 

Y e = -3C/ 13 mX 2 y 2 A 2 + y x 

\ y 9 X 2 -y 2 / V & / 
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^0 





2/6 A \ 


y 4 A 2 


£2 





A 2 + /3| 










^ 9 A 2 


-J/2 






Vo 








The {7 matrix defines the light SU(2)l Higgses in term of the PS Higgs field components, 
as explicitly written in appendix 



U21U 



(i) 



B 



and yi has to be read as {/ujff + Ui 2 y^ while as 



With respect to the mass matrices given in eqs. (47-49) we have reabsorbed a power 
of A. 



5.2 Analytical approximations 

In the appendix [E], we report all the RGEs for the Yukawa matrices, while here we discuss 
the results. The RGEs present the general compact expressions 



dt' 
dt' 



Tr. 



Yf>Y}, 



YfYj, 



Y L + Y L T T L 



YfY}, 



Y f + 



G f ^(Y r Y},)\-J249^ 



Y L 



Yt 



(86) 



where the index / runs over {e,u,d}, the parameter t' is defined as t! = t/(167r 2 ) = 
log/i/(167r 2 ), Tx [• • •] is a matrix written in terms of the fermion Yukawa matrices Yy, 
Gx [■ • •] is function of the trace in the flavour space over the Yukawa matrices Yf> and c{ 
are the Casimir coefficients related to the group representations (see appendix [E] for the 
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details). The generic solutions are given by 



yl(^) ~ Y\. e cViTi x exp 

i 

x Y L (fio) x exp 
Yf(ji) = Y\e~°i Xi x exp 
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t'(po) 
*'Oo) 



Tr(Y>r;,)" 














dt' 


x exp 


/ ^ 


d£ 


X 













XT 



dt' 



Tr(r f r/,)" 








Y f ,Y}, 






dt' 


x exp 


/ ^ 


dt' 


X ^(/Xo 






Jt'(po) 







(87) 

where Xj = J^^l gi(t') 2 dt'. When we fix /i ~ Mr and fi ~ Mst/sy these formulas can be 
approximated by 



Y L {M { 



SUSY; 



I + Gl 



Tr (Y f ,Y},)\ At' - £ c i x i Y L {M T )+ 



Yf'Yl Y L {M T ) + Y L {M T )Fl 



Y f (M : 



SUSY) 



f 

1 + 0, Tr (Yf/Yj, 



At' 



At'-J2^)Y f (M T )+T f 



Yf>Y}, 



Y f (M T ) At' 



with At' = l/l67r 2 \og(M SUSY /M T ) « 0.13. 

For what concerns the quark sector we find the following approximated expressions for 
the masses of the last two families and the Cabibbo angle 



m 2 t (M S usY) 


~ 


m 2 c (M SU SY) 




m 2 b (MsusY) 


~ m 6 2 (l 


m 2 s (M SU SY) 


~ m 2 a (l 


q 12 (Msusy) 


~ 0?2 + 



S9) 



(l-2E i cfX, + 16 7 Af) , 
(l-2£ iC ?X i + 8 7 A<) , 

J-U 2 3 XAt' , 

where 7 = m 2 /(v± + fiv^) 2 and the masses and the Cabibbo angle on the right of the 
previous expressions are intended at the Mt scale. Note that the demand that ml is still 
positive at Msusy gives an upper bound on 7 of 0.7. The charged lepton masses are very 
similar to the down-quark masses and indeed we have 



m 2 T {M S usY) 



ml(M : 



SUSY) 



™>l (1 - 2 EiCf2i + 16 7 At') 
m / 2 t (l-2£^X i + 8 7 At') . 



(90) 



We now consider the neutrino sector (see [36] for a general approach at RGEs with or 
without flavour symmetries) and the modification due to the RG running. We remember 
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that our model at the GUT scale naturally predicts the QD spectrum, with both normal 
and inverse ordering, while choosing a less natural range of the parameters space we may 
have both NH and IH spectrum. To analyze the effect of the RGEs on the neutrino mass 
matrix M v , it is worth rotate M u of a maximal rotation in the (12) sector. Then, at the 
high scale, M v is diagonal and reads as 

/ k' - Xk[ \ 

M v = J k' Q + Xk[ \v L , (91) 
\ k' + X 2 k' 2 J 

where without loss of generality k' Q can be taken real, by a redefinition on the phases. After 



the running at M S usy eq. (91) gets a correction AM V that is given by 



AM V 

vl 



Xk[ 






(92) 

/ -2A/c^ - 13A^7/2 -^7/2 + k[X^/2V2 \ 

+ At' -&o7/2 2Afc£, + 13Xk[j/2 k' w + + k[Xj/2V2 

\ k' w^ + k[Xj/2V2 k' w + + k[Xj/2V2 7~/k' ) 

with w± = (±7/y/2 + V2X), for the QD and IH case when k' 2 ~ 0(1). For the NH case 
characterized by k' 2 ~ A -2 , AM V assumes the following form 



/ 13 \ Xk 'i 

~K J2 c i Xi + T^ 7 A*' 1 + E dfli -Xk[ 
\ i A ) i \ 





1 ! 

-2A^ - 13Afci7/2 -^7/2 +k'_"//2V2 \ 

+ At' [ -&07/2 2A^ + 13AA; , 1 7/2 k' w+ + k' + j/2^/2 

k' w„ + k'_-f /2V2 k' w + + k' + j/2V2 7j(k^ + k' 2 X 2 ) J 



(93) 



with k' ± = (k' 1 X±k' 2 X 2 ). 

We can now consider the three different cases QD, NH and IH, which the model accounts 
for at the GUT scale. 

- QD case k' , k[, k' 2 ~ 0(1). 
The correction given by the running induces a rotation in the (23) sector characterized 
by 

tan2^ 3 ~ -^^7 ~ -2V2A^7, (94) 

being At' ~ 2A 2 . This is a large contribution, which deviates the atmospheric angle 
from the initial maximal value, spoiling the agreement with the experimental data at 
3a. A possible way-out would be if this large correction is erased by a corresponding 
large correction in the charged lepton mass matrix. However this is not the case, 
because the maximal 8 23 in the charged lepton mixing matrix is stable under the RG 
running. As a result, the QD case is not viable. 
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- NH case =► k' ~ A 2 , k[ ~ 0(1) and fc£ ~ A" 2 . 

The corrections both for the atmospheric and the reactor angles are of order At ,r y/2\/2 ~ 
3A 3 and can be safely neglected. Analogously, also the mass splittings receive devia- 
tions which can be neglected. On the other hand, the charged lepton mixing matrix 
is stable under the RG running. As a result the three mixing angles at M S usy can 
be well approximated with their initial values at Mt- 

- IH case k' Q ~ A 2 and k[, k' 2 ~ 0(1). 

All the corrections to the neutrino mixing are of order 7 At'/ \/2. While the solar mass 
splitting receives negligible contributions, the atmospheric mass splitting is corrected 
as follows 

Am 2 atm {M SUSY ) = Ami ^ " 2 E <^ + 13 7 At ') " ( 95 ) 

Combining the neutrino mixing with the charged lepton one we get the following 
mixing angles for the lepton sector: 

0[ 2 (M SUSY ) ~ tt/4 - 9\ 2 + ^ 3 7/2 At' , 

9 l 23 (M SUSY )~ ^3 + 7At', (96) 
0[ 3 (M SU sy)~ 6 e 13 -6f 2l /2 At', 

where the corrections coming from the RG running are proportional to At'. In par- 
ticular the corrections to 9[ 2 (Msusy) an d 9[ 3 (M S usy) come from the charged lepton 
sector, while that one to 9 1 23 (Msusy) arises only by the neutrino sector. This distinc- 
tion strictly holds only at this order of approximation, while considering contributions 
of the order of A 3 all the angles receive corrections from both the sectors. 

- IH case k 2 ~ A 2 , k' ~ 1 and k[ ~ A -1 . 

With respect to the previous case, the solar mass splitting gets a non-negligible 
correction, that can be written as 

Aml ol {M SUSY ) = Am 2 sol (1 - 2J2 i c^X i + 14 7 At') . (97) 

For what concerns the atmospheric mass splitting and the lepton mixing angles we 
recover the same shifts as in the previous case. 



6 Neutrino Phenomenological Analysis 

In the previous section, it was concluded that the quasi-degenerate spectrum is unstable 
under the RG running and becomes phenomenologically inviable due to too large cor- 
rections to the atmospheric mixing angle. Only the normal and inverse hierarchies were 



shown to be stable and phenomenologically viable. This is different from 23 , where only 



the normal hierarchy and the quasi-degenerate spectrum with normal ordering were found. 
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In this section, we will discuss neutrino phenomenology in more detail, focussing on the 
value of the reactor mixing angle 8[ 3 and the possibility of neutrinoless double beta decay. 



The analytical expression of the reactor mixing angle is given by eq. (64) at NLO and 



it is corrected by the RG running as in eq. (96) for the two IH scenarios studied in the 
previous section. We see that 9[ 3 is typically of order A ~ 0.2, so sin 2 9[ 3 ~ 0.04, which is 
rather large, but still allowed at the 3a level as shown in table [T} Here, we complete the 
study on 8[ 3 , performing a numerical analysis. As can be seen from eqs. (47)-(54), the 



neutrino and the charged lepton mass matrices at NLO are a function of many parameters. 
However the GUT nature of the model allows us to fix most parameters that occur at LO 
because they enter in the low energy expressions for quark and charged lepton masses, as 



can been seen comparing eq. (28) with eq. (89). Note that all dimensionless parameters, 
i.e. the can be fixed of order 1. 

The other free parameters can be fixed as random numbers of order 1, except for the 
cases where we have argued in the previous section that they should have slightly larger 
or smaller values. Because of the use of random numbers, the predictions of our model are 
no longer single valued. 
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Figure 3: The solar angle versus the reactor angle. On the upper line the two IH cases of the previous 
section (on the left the first case and on the right the second one), while on the lower line the NH one. 
The two vertical lines are the 3a bounds for sin 2 0[ 2 according to [l] . The upper horizontal line is the 3a 
upper bound for sin 2 6[ 3 and the middle line is the best fit value. The lower line is the 95 % exclusion 
confidence level after 3 years of Daya Bay data taking (5). 
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IS 



We plot the reactor angle versus the solar angle in figure |3j At NLO, eq.(65), 9[ 2 
driven away from the maximal value 7r/4 by a term proportional to A (note that we take 
only the corrections which decrease the value of the solar angle, neglecting those which 
increase it). We see that this deviation is not for all values of the parameters large enough 
to bring it in the observed region, although this happens for a significant number of them. 
As explained above, larger values of sin 2 9[ 3 are favoured and almost all points are in the 
sensitive region for experiments. 

To study neutrinoless double beta decay, we consider the effective Ou/3/3 parameter m ee , 
defined as 

m ee = [U diag(mi, m 2 , m 3 ) U] n . (98) 

In figure [4] we plot m ee against the lightest neutrino mass, which is mi and 777,3 i n the NH 
and IH case respectively. The future experiments are expected to reach good sensitivities: 



90 meV (38J (GERDA), 20 meV (39J (Majorana), 50 meV [40] (SuperNEMO), 15 meV (41 
(CUORE) and 24 meV (42) (EXO). As a result, looking at figure [§ the whole IH band will 
be tested in the next future and with it the two cases of our model which allow for the IH 
spectrum. 
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Figure 4: Neutrinoless double beta decay plots. On the upper line the two IH cases of the previous section 
(on the left the first case and on the right the second one), while on the lower line the NH one. The 
background red (blue) points refer to the allowed region for the NH (IH), taking into account the lepton 
mixing angle values with their 3(7 errors. 
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7 Conclusion 



In this paper we have addressed several aspects of the interplay between a GUT based 
model and a discrete flavour symmetry. The paper should indeed be considered as the 
combination of two distinct parts: in the first one we mainly discussed the building of the 
model from the flavour point of view, while in the second one we faced the problem to 
justify the assumptions made in the first part and to achieve the correct gauge symmetry 
breaking chain. 

More in detail, the symmetry group of our model is PS x Gf, where PS stands for the 
GUT Pati-Salam gauge group SU(4) C x SU{2) L x SU{2) R and Gf for the flavour group 
5*4 x Z4 x U(1)fn x U(1)r. Within this GUT context one has the relationship between 
the down-quark and charged leptons mass matrices, Md ~ M e , which can easily be used 
to revise the old idea of quark-lepton complementarity. In the model this is obtained 
by the use of the non-Abelian discrete flavour symmetry S4 properly broken through the 
VEVs of a set of flavon fields, which transform as triplets under 5*4. The additional Abelian 
symmetries, which enters in G/, play different roles: Z 4 keeps quarks separated from leptons 
and neutrinos from charged leptons and prevents dangerous couplings in the superpotential 
of the model; U{1)fn helps to justify the charged fermion mass hierarchies; U{1)r is a 
common ingredient of supersymmetric flavour models. It contains the discrete .R-parity 
and is useful to build a suitable flavon superpotential that allows the correct S4 breaking 
pattern. 

Already at the leading order, the model shows nice features: we are able to reproduce 
the mass hierarchy between the third and the second charged fermion families, the bottom- 
tau unification, the Georgi-Jarlskog [30] relation \m^\ = 3\m s \ and, under the assumption 
of type-II See-Saw dominance at the GUT scale, a realistic neutrino spectrum. However at 
this level of approximation, both the CKM and the PMNS mixing matrices are not correct: 
the quark mixing matrix coincides with the identity matrix, while the lepton one is given 
by the BM pattern. It is worth to recall here that the BM mixing corresponds to maximal 
solar and atmospheric angles and to a vanishing reactor angle: only the solar angle is not 
in agreement with the data as it deviates from the experimental central value by a quantity 
close to the Cabibbo angle, A ~ 0.2. 

At next-to-leading order, the wrong predictions for the fermion mixing angles are cor- 
rected: in the CKM matrix, the mixing angle 9\ 2 receives contributions of the order of A, 
fitting the value of the Cabibbo angle; analogously, in the PMNS matrix, the solar angle is 
corrected by the same amount and we find the nice result that 8[ 2 ~ 7r/4 — O(X). At the 
same time, also the reactor angle receives significant contributions and indeed at this level 
of approximation it results 9[ 3 ~ 0(A): this is an interesting feature of our model, because 
this value is close to its present upper bound and it will be tested in the forthcoming 
neutrino experiments [4jj9]. 

Once we consider the higher order terms, we find the other two CKM angles of the 
correct order of magnitude, 6 23 ~ 0{\ 2 ) and 6f 3 ~ 0(X 3 ), and small corrections are 
introduced in the PMNS angles: in particular the atmospheric angle becomes 8 l 23 ~ 7r/4 + 
0(X 2 ), justifying a small deviation from the maximality. For what concerns the masses, 
all the fermions are massive and the mass hierarchies fit the experimental observations. 

On the other hand, the neutrino spectrum could be either quasi degenerate or normal 
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or inverse hierarchical. Only the first case corresponds to a completely natural choice of the 
parameters, which, in the absence of an explanation coming form a higher energy theory, 
should be of order 1: in order to allow the NH and the IH, the parameters should span in 
a larger range of values, namely A~ 2 -j- A 2 . 

In the second part of the paper we have studied the Higgs scalar potential and the 
running of both the gauge couplings and the Yukawa mass matrices under the RGEs. With 
this analysis we looked for the constraints which arise to justify the Higgs field VEV pattern 
used in the flavour section and the assumption of the type-II See-Saw dominance. The 
presence of the flavour group Gf modifies the Higgs field content necessary to implement 
the classical breaking pattern of the PS gauge group and as consequence not all the results 
obtained by studying minimal versions of PS are recovered. In particular we need the 
presence of two PS multiplets (15, 1, 1), A and B, responsible to break the unified colour 
symmetry SU(4)c to SU(3)c x U{1)b-l- Two copies of the (1,2,2) multiplet, and 0', 
and one (15, 2, 2) field, p, are necessary to implement the condition v u p = . Lastly, we need 
the new fields E, £' ~ (1, 3, 3) and £ ~ (1, 1, 1) to have a type-II See-Saw contribution at 
tree level. Gauge coupling runnings are affected by the large field content and in particular 
we found that the requirement of having type-II See-Saw dominance constrains the model 
in a small region of the parameter space, in which all the heavy mass scales are sandwiched 
between 10 11 GeV and 10 13 GeV. At the same time Yukawa mass matrix running shows 
that while the CKM Cabibbo angle is stable under the RGEs evolution, the PMNS mixing 
angles are stable only if neutrinos present a NH or an IH spectrum, ruling out the QD 
case. As already stated, the QD spectrum would be the most natural and probable case at 
the GUT scale, but the Yukawa RGEs analysis further reduces the allowed region of the 
parameter space. 

In section [6j we performed a brief phenomenological analysis on neutrino observables 
considering all the constraints which come from the flavour and the Higgs sectors. We have 
first considered the value of the reactor angle in terms of the deviations of the solar angle 
from the maximal value and the numerical analysis confirmed the analytical results: in 
our model, 6[ 3 naturally acquires a value not far from its present upper bound. After this, 
we studied the neutrinoless double beta decay effective mass, m ee , and we have seen that 
the next future experiments are expected to reach sufficiently good sensitivities to test our 
model in the IH regime. 

As a final comment, we underline that there is a strong tension in combining a GUT 
model with a (discrete non-Abelian) flavour symmetry and therefore a parallel study is 
not only interesting but also recommended to provide a viable model. In the first part 
of the paper we produced a flavour-GUT model that is quite realistic and viable. Only 
the study done in the second part, regarding the Higgs sector reveals that the model is 
restricted into a small region of the parameter space, reducing the freedom in the choice 
of the parameter values. Even if our results are model dependent, our construction shares 
many features with other models present in literature, where often a detailed discussion 
of the Higgs sector is missing. In our opinion, this is a serious drawback and we would 
suggest to consider the interplay between GUTs and flavour symmetries in this kind of 
models as well. 
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A The symmetric group S4 



In this appendix we report the character table and the Clebsch-Gordan coefficients of the 
S4 discrete group in our basis. 





n 


h 


Xii 


Xu 


X2 


X3i 


X3 2 


Elements 




1 


1 


1 


1 


2 


3 


3 


1 


c 2 


3 


2 


1 


1 


2 


-1 


-1 


T 2 , ST 2 S, ST 2 ST 2 


c 3 


6 


2 


1 


-1 





1 


-1 


S, T 3 ST, TST 3 , T 2 ST 2 , ST 2 ST, TST 2 S 


c 4 


8 


3 


1 


1 


-1 








TS, ST, (TS) 2 , (ST) 2 , T 2 ST, TST 2 , T 3 ST 2 , T 2 ST 3 
T, T 3 , ST 2 , T 2 S, STS, TST 


c 5 


6 


4 


1 


-1 





-1 


1 



Table 7: Character table of the S4 discrete group, d are the conjugacy classes, n the number of elements 
in each class, h the smallest value for which x h — 1 • In the last column we have reported the elements for 
each class in terms of the generators of the group. 

The generators, S and T, obey to the following rules 

T 4 = S 2 = (ST) 3 = (TS) 3 = 1 , (99) 
and are of the following form for the five different representations: 
li: 5 = 1, T 
1 2 : S = -l, T 

2 : S = '- \ 'r, v ; I . T 




3i : S = 1 , T 



5=0-10, T 




Using these generators we calculate the Clebsch Gordan coefficient for all the Kronecker 
products. In the following we use ctj to indicate the elements of the first representation of 
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the product and fti to indicate those of the second representation. 

We start with all the multiplication rules which include the 1-dimensional representations: 

11 <8> r] = rj ®1\ = rj with i] any representation 

12 <8> I2 = li ~ a/3 



2=2 



3i 



3i 



-a/3 2 

a/3 2 
a/3 3 

a(3 2 
af3 3 



The multiplication rules with the 2-dimensional representation are the following: 

li ~ etipi + a 2 (3 2 



2©2 = liffil 2 ©2 with < 



1 2 ~ -aip2 + 
2 - 



2®3i =3i©3 2 



with < 



2 <g> 3 2 = 3i © 3 2 



with < 



«i/3 2 + a 2 /3i 
otifii - a 2 (3 2 



_a 2 /?i 

3 j I -|(-\/3q;i^2 + tt2/?2) 
l(y/3aip 3 - a 2 /3 3 ) 

|(a/3q;2^2 - (X1P2) 

-\{\fZa 2 f3 3 + ai/33) 
-|(-n/3Q!2^3 + en/33) 



\(\ / 3a 1 f3 3 - a 2 (3 3 ) 
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The multiplication rules with the 3-dimensional representations are the following: 



3i ® 3i = 3 2 ® 3 2 = li © 2 © 3i © 3 2 



3i 



3i 



li ~ ai/5i + a 2 (3 2 + a 3 (3 3 

-^(a 2 (3 2 - a 3 (3 3 ) 



with 



with 



3i 



-2ai^i + a 2 /3 2 + a 3 (3 3 ) 



3i 



ai/fe + a 3/3i 
«i/3 2 + a 2 /3i 

a 3 /3 2 - a 2 /3 3 
«i/3 3 - a 3 /3i 
a 2/ 5i - ai/3 2 



1 2 ~ + a 2 /3 2 + a 3 /3 3 

^(2ai/3i - a 2 /3 2 - a 3 /3 3 ) 
^(a 2 (3 2 - a 3 /3 3 ) 



"3^2 - "2/^3 

ai/33 - a 3 /3i 
a 2 /3i - «i/3 2 

a 2 /?3 + a 3 /?2 
ai/?3 + a 3 /?i 
ai/3 2 + a 2 /3i 



B Higgs scalar spectrum 

We now proceed to presenting the scalar mass matrices for all the fields introduced in sec. [4] 
according to the SM~ SU(3)c x SU{2)l x U(1)y representations, indicating their origin 
with respect to the PS and the colour-broken Pati-Salam (CbPS) phase. We recall that in 
the colour broken phase the symmetry group is given by SU(3)c x SU{2) L x SU{2) R x 
U(l) b-l- For completeness for each field we also indicate the corresponding T 3R value, with 



T 3 r the diagonal generator of SU(2)r. We use the same notation as in 37 , in which we 



write the Dirac scalar mass matrices as they could be read directly from the superpotential 



of eq. (76) at the scale Mr. We label the mass matrices with S, D and T according to the 



singlet, doublet or triplet representations of the SU(2) L gauge symmetry. 
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B.l SU(2) L singlets 
B.l.l Singlets (1, 1,0) 



/ 



M, 



si 













^X R M R 



f 2 \ R M R f 2 \ R M R 



V2 

f 2 \ R M R 

















V2M Cl X A + x X AB ^ + j^X B M C2 j^X ABi M C2 

^XbzM Ci - y ^Xab(M C2 
^X AB ^M Cl 



\ 



X AB ^ + ^X B Mc 2 



M c ,Mc 2 x 

(100) 



with 



x = 



A\ A M 2 Ci - 2V2X ABi M C2 V i - 2X b Mq 2 - 8X R M 2 R ) 



2V2M Cl 

y = -^-(^2\ AB ^ + 2X B M C2 
V 2M C2 



(101) 



Msi has a vanishing eigenvalue that is eaten by the corresponding gauge boson. Moreover 
it can be checked that one of the singlets, which we call £o, has a mass ~ M R while all the 
others masses appear as combination of Mqi, M C 2 and V^. 
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B.l. 2 Singlets (3, 1, 2/3) © (3, 1, -2/3) 
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I V2X A M Cl + x Xab^ + -^X B M C2 

X AB ^ + ^X B M C2 ^X B M Cl - y 
\ 2X R M R 



2X R M R \ 


-V2X R M Cl J 



(102) 
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M S 2 has a vanishing eigenvalue: it corresponds to the massless GBs (3, 1, 2/3) ©(3, 1, —2/3) 
eaten by the gauge bosons. 
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B.1.3 Singlets (8,1,0) 



1 V(\abc-^\bM C 2 -^MciXb-V 
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B.1.4 Singlets (1,1, ±1) 



These states correspond to the last two massless GBs eaten by the respective gauge bosons. 
Indeed with these last two massless states the total amount of GBs eaten is 9 = 1 + 2x3 + 2, 
as it should be for the breaking SU(A) C x SU(2) L x SU(2) R ->■ SU(3) C x SU(2) L x U(1) Y - 
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B.1.5 Singlets (1,1, ±2) 

Two massless charged singlet (5 ++ and 5 ++ ). 
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B.1.6 Singlets (3, 1, -1/3) © (3, 1, 1/3) 

-V2M C1 \ R 
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B.1.7 Singlets (3, 1, -4/3) © (3, 1, 4/3) 

-V2M C1 X R 
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B.1.8 Singlets (6, 1, 4/3) © (6, 1, -4/3) 

-2V2M C1 \ R 
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B.1.9 Singlets (6, 1, 1/3) © (6, 1 - 1/3) 

-2V2M C1 X R 







PS 


CbPS 






A* 


(10,1,3) 


(6,1,3,2/3) 





i?l 


A* 


(TO, 1,3) 


(6,1,3,-2/3) 






B.1.10 Singlets (6, 1, -2/3) © (6, 1, 2/3) 

-2V2M C1 X R 
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B.2 SU(2) L doublets 
B.2.1 Doublets (1,2, ±1/2) 



M, 



Dl 



M H , My v ; ,.U r2 A,, v 



When M| )1 is diagonalized according to 



U T ■ M 2 m ■ U = M m M, 



(104) 



DUV1D1 



(105) 
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we get three up-type (down-type) Higgs doublets h u , h' u , H u (hd, h' d , Hd). 

(106) 



/ 4>u,d \ 






4>' u , d 






\ Pu,d J 




\ H u,d ) 



Therefore for the up (down) projections we have 



u.d i± \ tt u,d , j j u.d , j j u.d 

V = {4> u ,d) = U n vi + U 21 v 2 ' +U 31 v 3 ' , 

v U / = Kd) = U 12 v u { d + U 22 v u / + U 32 v^\ (107) 

v u / = ( Pu ,d) = W + W + U 33 v u / , 

where = (h U7 d) -,v^ d = (h' ud ) ,v^ ,d = (H u ^)- In general a light doublet, i.e. massless 
at the Mc scale, gets a mass term at Msusy and its VEV, at the EW scale, is of order of 
the EW scale ~ vw- On the contrary for a heavy doublet of mass M, its induced VEV 
at the EW scale is ~ v^y/M, that for M ~ Mc is completely negligible with respect to 
vw Consider now the condition (p u ) = that we imposed to get the correct fermion mass 
matrices. In the standard case we would have only one up- and one down-type light Higgs 



doublets, being all the other doublets heavy. Assume now that h u d in eq. (106) are the up- 



and down-type light doublets at Mc- Then we have v ~ vw while v%' d ~ v%' ~ 0. From 



eq. (107) we see that in this case it would be impossible to make vanishing the p projection 
along the up-direction (that implies U\ 3 = 0) still maintaining a non-vanishing (pd) p]For 
this reason the condition (p u ) = implies a non-standard scenario and the presence of 
two light doublets of up-type at the Mc scale, namely h u and h' u . However the symmetric 
nature of Mdi ensures that as consequence we are also left with two down-type light Higgs 
doublets, hd and h' d . Nevertheless this does not imply that (pd) vanishes because we have 

(Pu) = U 13 v^ + U 23 v^ = 0, 

(Pd) = M + M^> (108) 

since v^ ,d depend on the soft terms. 

In conclusion we have to impose two constraints on the free parameters that enter in 
Mdi corresponding to require that Mdi has two vanishing eigenvalues. Notice that the 
condition of having Md 1 of rank 1 is fine-tuned but it is not more fine-tuned than imposing 
Mdj of rank 2, which is universally accepted whenever the MSSM has to be recovered. In 
our case the fermion mass matrix structures impose a slightly different condition —M^ of 
rank 1 - but both the requirements are satisfied by fine-tuning the parameters that enter 
in the mass matrix. 



u It is worth to say that even the condition [7 13 = is not natural and easily realized. Therefore in the 
most general case with only one up-type (down-type) light doublet, the condition (p u ) = implies that all 



the other VEVs given in eq. (107) vanish 
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B.2.2 Doublets (3, 2, 7/6) © (3, 2, -7/6) 
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B.2.3 Doublets (3, 2, 1/6) © (3, 2, -1/6) 
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B.2.4 Doublets (8, 2, 1/2) © (8, 2, -1/2) 

M p - -^=\p A M Cl 
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B.3 SU(2) L triplets 

B.3.1 Triplets (1, 3, 1) © (1, 3, -1) 

/ M s \ ( Vz \ 

M T1 = hV ( M s , (109) 
V j= 2 XM R M AL + ^X L M Cl J 

with 
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B.3.2 Triplets (1,3,0) 
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B.3.3 Triplets (3, 3, -1/3) © (3, 3, 1/3) 
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B.3.4 Triplets (6, 3, 1/3) © (6, 3, -1/3) 



M Al - -^\ L M C i 



with 







PS 


CbPS 


T3R 




A L 


(TO, 3,1) 


(6,3,1,-2/3) 





Ri 


A L 


(10,3,1) 


(6,3,1,2/3) 






C NLO contributions to the flavon scalar potential 

The superpotential Wd, linear in the driving fields Dr, <^r, xr and <jr, is modified into: 

w d = w° d + 5w d . (Ill) 

where 5wd contains the NLO contributions, suppressed by one power of 1/A with respect 
to Wd- The corrective term Swd is given by the most general quartic, S4 x ^-invariant 
polynomial linear in the driving fields, and can be obtained by inserting an additional 
flavon field in all the LO terms. The Z 4 -charges prevent any addition of the flavons (p and 
<// at NLO, while a factor of a or x can be added to all the LO terms. The full expression 
of 5wd is the following: 



6wd = i (E x ^ R + E w ^ R + E s ^ R + E v ^ R ) 

\i=l i=l i=l i=l / 



(112) 



where Xi, Wi, Sj and Vi are coefficients and /j Dii , represent a basis of 

independent quartic invariants: 



I° R = GrGGO 

I a 2 R = aMxxhi) 



= w(xx) 



(113) 
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In the previous terms we indicate with (. . .) the singlet li, with (. . .)' the singlet I2 and 
with (. . .) R (R = 2, 3i, 3 2 ) the representation R. 

The NLO flavon VEVs are obtained by imposing the vanishing of the first derivative 
of Wd + Swd with respect to the driving fields a R , \Ri D R and ip R . We look for a solution 
that perturbs eqs 
the flavons $ 



(17) and (18) to first order in the 1/A expansion: for all components of 
we denote the shifted VEVs by 



($) = ($) LO + 5$ 



(117) 



where ($)lo are given by eqs. (17) and (18). 



It is straightforward to verify the following results. In the Majorana mass sector the 
shifts 5a, 5\ turn out to be proportional to the LO VEVs ($)lo and can be absorbed in 
a redefinition of the parameters v x and v a . Instead, in the Dirac mass sector, the shifts 
Sip, Sip' have a non-trivial structure, so that the LO texture is modified: 



(<P) 



Sv u 



<p 
,1 



Sv 



:n8) 



-v, 



where v' and v' , satisfy a relation similar to that in eq. (73) and the shifts Sv^, and Sv 



V 



are suppressed by a factor A with respect to the LO entries v' and vL, respectively. 



D Gauge coupling running 

In this appendix, we provide the coefficients of the /3-functions for the gauge coupling 
running in the different regimes. The complete matter fields run from the GUT scale 
down to the Msusy scale, where the SUSY partners decouple. For what concerns the 
scalar fields in section [4], we have already outlined the scalar spectrum according to the 
different scale at which the fields decouple. As a result the computation of the /3-functions 
is straightforward. Calling \x the generic scale, we have 

- for Mc < fi < Mqut'- all matter is in the left and right handed multiplets (4, 2, 1) 
and (4, 1,2) as mentioned in table [2] In the Higgs sector, we have all the fields 
mentioned in table [U This leads to the coefficients 

Psuwc = 54 , SU (2) L = 69 , (3 SU(2)R = 69 . (119) 
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Due to the large matter content these coefficients are very large and the /3-functions 
are very steep. As consequence the theory is in the Pati-Salam regime only for a very 
small range of energies, as can indeed be seen in figure [2j Almost after passing the 
scale M c , the SU(2) R coupling constant enters the non-perturbative regime. 

- M R < fi < M c : the theory undergoes to the SU(3) C X SU(2) L x SU(2) R x U(1) B -l 
symmetry. Considering the matter (left and right handed doublets of quarks and 
leptons characterized by different U(1)b-l charges) and the scalar fields, we find the 
coefficients 

/3su(3) c = -3 > Psu(2) L = 18 , Psu(2) R = 18 , f3u(i) B _ L = 24 . (120) 



- Mt < fj, < Mr. in this regime we have all the usual MSSM matter particles, four light 
higgs doublets (two up-type and two down-type), a couple of SM triplets (1,3, 1) © 
(1,3,-1) and two extra charged singlets (<5 ++ and 5 ++ ). The coefficients of the 
/3-functions are 

(3su(3) c = -3 , Psumz. = 4 , Pu(i) Y = 69/5 . (121) 

The hypercharge that appears in the last term is related to SU(2)r and the B — L 
charges in the previous regime by 



- Msusy < /i < Mt'- in this regime we have all the usual MSSM matter particles, four 
light higgs doublets (two up-type and two down-type) and two extra charged singlets 
(5 ++ and 5 ++ ). The /3-function coefficients are 

AsT/(3) c = ~3 , Psu(2) L = 2 , Pu(i) Y — 12 • (122) 
This should be compared with the (-3, 1, 33/5) coefficients of the ordinary MSSM. 

- vw < /i < Msusy'- we have the particle content of the standard model, with the 
exception that there are four Higgs doublets. We have therefore the following /?- 
function coefficients 

0su{3) c — ~7 > Psu(2) L = -8/3 , Pu(i) Y = 22/5 . (123) 
This should be compared with the (-7, -19/6, 41/10) coefficients of the ordinary SM. 



E Yukawa running 

As we already reported in section |5j because of the closeness of the intermediate scales 
between Mqut and My we can consider only the running between M? and Msusy to 
provide analytical approximations for the evolution of fermion masses and mixing under 
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the RGEs effect. At Mt the scalar SU(2)l triplet has been already integrated out giving 
rise to the effective Weinberg operator given by 



aijY Lr 



L r L s h Ui h U j 

Mr 



(124) 



with h Ul = h u , h U2 = h' u and coefficients arising by the scalar potential and Y L given in 



eq. (83) 



For what concerns the charged fermion Yukawa, at Mt the Dirac part of the superpo- 
tential is written as 

Y U Q U c h u + Y' U Q U c h' u + Y d Q D c h d + Y' d Q D c h' d + Y e L E c h d + Y' e L E c h' d , (125) 



with the Yukawa mass matrices given in eq. (85). 



The Yukawa matrices RGEs are therefore given by 
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dt 
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2 dY L 



dt 



Y' 

U 1 



3 Y u Yl + Y d Y\ + 3 YX ] + + 3 Tr(Y u Yj) - Z^gf Y„ 

3 Y u Yl + Y d Y\ + 3 YX ] + + 3 Tr^Y^) - ^cfg 

Y U YJ + 3 Y d Y\ + YX ] + 3 FJY^ + 3 Ti(Y d Y d ) + Tr(Y e Y £ t) - E <C fe 
Y w Yj + 3 Y d Yi + Y^t + 3 YjY^ + 3 Tr (YjY^) + Tr(Y e 'Y e 't) - E iC ^ 2 



Y' 



3 Y e Y e t + 3 Y e 'Y e 't + 3 Tr (Y d Yj) + Tr(Y e Y e t) - £,c^ 2 ] Y e 
3 Y e Y e t + 3 YX ] + 3 Tr(Yp?) + T^Y^) - ^Adl 



YL. 



[(Y e Y e t + fx^l + n(nn f + n%T - ^gf)y L ] («u + « 12 + a 22 ) 

+ [6 Tr(Y u Yj> n + 3 Tr(Y u Y M t)a 12 + 3 Tr(Y^Y^)ai2 + 6 Tr^Y^)^] Y L , 

(126) 
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